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(54) Method for determining a partition in order to insert a watermark, and associated insertion 
and decoding methods 



(57) A method of determining a partition into distinct 
regions of a set of coefficients representative of at least 
a part of an image in order to insert a watermark com- 
prises the following steps: 

generation (E501 ) of a centred pseudo-random se- 
quence (w) equal in size to the cardinal number (N) 
of the said set of coefficients, formed from centred 
pseudo-random sub-sequences; 
modulation (E502) of the said set of coefficients by 
the said centred pseudo-random sequence (w) in 
order to insert the same information bit on the said 
set of coefficients; and 

checking (E504-E508), over each region, of a crite- 
rion of detectability of the information bit inserted by 
modulation. 

Use especially for determining an adaptative parti- 
tion as a function of a criterion of detectability of an in- 
formation bit inserted on each region. 
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Description 

[0001] The present invention concerns a method of determining a partitioning of a digital image tor inserting a wa- 
termarking signal, a method of inserting a watermark, a method of processing an image and an associated method of 
decoding a watermarking signal. 

[0002] Correlatively, it concerns a device for determining a partitioning of a digital image for inserting a watermarking 
signal, a device for inserting a watermark, a device for processing an image and a device for decoding a watermarking 
signal. 

[0003] The present invention lies in general terms in the technical field of watermarking of digital images and more 
particularly still images. 

[0004] Watermarking digital data makes it possible to protect these data, for example by associating copyright infor- 
mation with them. 

[0005] In its general principe, watermarking consists of inserting an indelible watermark in digital data, which can be 
likened to encoding additional information in the data. 

[0006] Decoding this additional information makes it possible to check the additional information which has been 
inserted. 

[0007] This inserted watermark must consequently be both imperceptible, robust to certain distortions applied to the 
digital image and reliable to detect. 

[0008] Conventionally, a usual technique of inserting a watermarking signal in a digital image consists of using a 
linear modulation model in which at least one subset of coefficients representing the digital image is modulated ac- 
cording to this linear model using a weighting coefficient. 

[0009] By denoting as X = {X i( 1 < i < N} a set of the coefficients representative at least of part of a digital image, and 
as w = {w jf 1 < j < P} a watermark of size P < N, a pseudo-random signal with known distribution and average of zero, 
the linear insertion formula is : 

XV = Xj + b ccj w., with 1 < j < P, 

in which {X jf 1 < j < P} is a subset of the set of coefficients X, b is an information bit, and a is a weighting coefficient 
30 also called modulation amplitude. 

[0010] The detection of the watermark then consists in detecting whether or not the pseudo-random sequence w 
has been inserted into a set of coefficients. This detection is carried out without using the original image and can be 
based on a standardised statistical test which makes it possible to calculate a probability of correct detection. 
[0011] Such an insertion technique makes it possible, by insertion of a watermark, to insert a single information bit 
35 since the response from the detector is binary (yes/no). 

[0012] In order to insert a larger number of information bits into the digital image, in particular when a code of Q bits 
is wanted, indicating, for example, the name or the address of the owner or of the author of the image, it is necessary 
to reiterate the previously described insertion method as many times as there are information bits to be inserted. 
Typically, in order to insert a binary signal, either b = 1 , or b = -1 is used. 

[0013] Put another way, Q subsets of coefficients have to be chosen and the modulation of these subsets has to be 
carried out by choosing Q watermarks. 

[001 4] Separate subsets of coefficients are preferably chosen such that the modulation operations are not superim- 
posed on one another, which could disturb the detection or cause troublesome visual effects. 

[0015] It is a matter, consequently, of choosing a partition of the coefficients representative of the digital image into 
Q separate subsets, each carrying one information bit. 

[0016] Numerous known methods use a technique of inserting a watermarking signal of given size by spectrum 
spreading. The disadvantage of these methods, describe for example in the article entitled "Secure spread spectrum 
watermarking for multimedia" by I.J. COX et al, in Proc. ICJP, pages 243-246, Sept. 1996 and in the article entitled 
"Digital watermarking of raw and compressed video" by F. HARTUNG et al, in Proc. SPIE 2952 : Digital Compression 
Technologies and Systems for Video Communication, pages 205-213, October 1996, is that they use an arbitrary 
partitioning of the image into blocks of fixed size without any guarantee of detectability of the modulations effected on 
each of the blocks. 

[0017] The aim of the present invention is to propose a method of determining a partitioning of a digital image which 
is adapted to the image, ensuring detectability of the bits inserted for a technique of insertion by modulation. 
[0018] The present invention thus relates to a method of determining a partitioning into distinct regions of a digital 
image in order to insert a watermarking signal by modulation of coefficients of said regions representing a digital image. 
[0019] In accordance with the invention, this partitioning into regions is effected by an adaptive partitioning as a 
function of a criterion of detectability of an information bit inserted in each region. 
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[0020] In this way an adaptive partitioning is performed which has the advantage of adapting to the content of the 
image signal and guaranteeing detectabiiity of the inserted watermarking signal, unlike the choice of an arbitrary par- 
titioning inio fixed blocks. 

[0021] This spatial adaptivity is necessary so as the content of the image is not spatially homogenous. 
[0022] This is because, for a given modulation amplitude, the greater the length of the sequence of the watermarking 
signal inserted, and consequently the greater the size of the region and therefore the number of modulatable coeffi- 
cients, the more reliable is the detection. 

[0023] It can in fact be shown that, when a correlation calculation is used for effecting the detection, the response 
of the detector is proportional to P 1/2 where P is the length of the sequence of the watermarking signal. 
[0024] Moreover, the response of the detector is also dependent on the image signal itself. The greater the variance 
in the coefficients representing the image, the weaker the response of the detector. 

[0025] In conclusion, for a given modulation amplitude and image, there exists a minimum length of the watermarking 
signal necessary for ensuring a given detection probability, this minimum length also depending on the local content 
of the image signal in the region in which the watermarking signal is inserted. 

[0026] The adaptive partitioning effected by the method according to the invention guarantees the detectabiiity of 
the watermarking signal inserted in the different regions determined by this partitioning. 

[0027] Preferably, the detectabiiity criterion is calculated using a maximum value of the weighting coefficient ensuring 
imperceptibility of the watermarking signal for the modulation of the set of coefficients representing said digital image. 
[0028] This is because, the greater the modulation amplitude, that is to say the value of the weighting coefficient, 
the more reliable is the detection in that a better probability of detection can be obtained. Nevertheless, the greater 
this modulation amplitude, the more this modulation may be visible in the watermarked image. 

[0029] It is therefore judicious to choose a modulation amplitude equal to a maximum amplitude, generally denoted 
JND (Just Noticeable Difference), beyond which the modulation becomes visibly perceptible in the image. 
[0030] It can be shown that this maximum amplitude also depends on the length of the sequence of the watermarking 
signal and that it decreases with this length. 

[0031] By choosing a maximum value of the weighting coefficient ensuring imperceptibility of the watermarking signal 
for the modulation of the set of coefficients representing said digital image, the value of the weighting coefficient on 
each region of the partitioning, and with a size less than the cardinal number of all the coefficients representing the 
image, will be much less than the maximum modulation amplitude which guarantees imperceptibility of the modulation 
on this region. 

[0032] The partitioning thus performed also ensures invisibility of the modulation effected on each of the regions 
obtained by the partitioning. 

[0033] According to an advantageous characteristic of the invention, the detectabiiity criterion is calculated using, 
for each coefficient to be modulated, a weighting coefficient according to a law of the form ctj = k^, where k t is a 
modulation factor dependent on the coefficients situated in a vicinity of the coefficient considered on said region and 
a v is equal to said maximum value of the weighting coefficient. 

[0034] It is thus possible to exploit the possibilities of visual "masking" of the modulation by the image signal itself. 
Typically, use is made of the fact that, the greater the activity of the image signal, in the sense for example of a significant 
local variance in the coefficients, the less visible will be the modulation. 

[0035] It is thus possible to locally increase the weighting coefficient and therefore the modulation amplitude, to the 
benefit of detectabiiity, without impairing the imperceptibility of the inserted modulation. 

[0036] According to another preferred characteristic of the invention, the method of determ ining a partitioning includes 
a step of applying a distortion before calculating the detectabiiity criterion. 

[0037] According to another preferred characteristic of the invention, the method of determining a partitioning includes 
a prior step of comparing the size of each region of the partitioning with a minimum size corresponding to the minimum 
size of a statistically significant sample for detecting a watermarking signal inserted in said region. 
[0038] Indeed, because of the calculations of a statistical nature performed when the watermarking signal is detected, 
notably when standardised hypothesis tests are used, it is necessary to consider a minimum boundary for the size of 
each region of the partitioning guaranteeing the validity of the detection tests. 

[0039] According to one particularly practical embodiment of the invention, for each region of the partitioning, all the 
coefficients of said region are modulated by inserting the watermarking signal, a detectabiiity amplitude is calculated 
from said detectabiiity criterion and said detectabiiity criterion is validated by comparing said detectabiiity amplitude 
with a predetermined threshold value. 

[0040] The detectabiiity amplitude can be the result of a calculation of correlation between the set of modulated 
coefficients and the watermarking signal or the result of a standardised statistical test such as the one described in 
the article entitled °A method for signature casting on digital images" by I. PITAS, in Proc. ICIP, pages 21 5-218, Sep- 
tember 1 996. 

[0041] Practically, the detectabiiity criterion can be validated if a minimum modulation length calculated for ensuring 
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detectability of the inserted watermarking signal in less than or equal to the cardinal number of said region 
[0042] According to a preferred embodiment of the invention, the method comprises the following steps: 

- generation of a centred pseudo-random sequence equal in size to the cardinal number of a set of coefficients 
representative of at least a part of the image, formed from centred pseudo-random sub-sequences; 

- modulation of the said set of coefficients by the said centred pseudo-random sequence in order to insert the same 
information bit on the said set of coefficients; and 

- checking, over each region, of a criterion of detectability of the information bit inserted by modulation. 

[0043] It is thus possible to achieve a partition of a set of coefficients into distinct regions, by carrying out an adaptive 
partitioning on the basis of a criterion of detectability of an information bit inserted on each region. 
[0044] Such an adaptive partitioning exhibits the advantage of adapting to the content of the image signal and of 
guaranteeing the detectability of the watermark inserted on the different regions determined by this partition, in contrast 
to the choice of a random partitioning. 

[0045] In practice, for each region of the partition, a criterion of detectability of an information bit inserted by modu- 
lation of the coefficients of the said region is checked. This criterion may consist in the comparison of a result of a 
standardised statistical test with a threshold value corresponding to a probability of detection fixed in advance. 
[0046] By virtue of the method in accordance with the invention, it is possible to generate a centred pseudo-random 
sequence on the basis of the merging of centred pseudo-random sequences of smaller sizes, and to modulate all the 
available coefficients in order to achieve the insertion of a watermark. 

[0047] By modulating all the coefficients with the same information bit and a centred pseudo-random sequence which 
is broken down into pseudo-random sub-sequences which are also centred on each region of the partition produced 
over the set of coefficients, the steps of generating a pseudo-random sequence and of modulation do not have to be 
reiterated at each new region tested. 

[0048] According to one preferred characteristic of the invention, the method of determining a partition comprises a 
single step of applying a distortion to the said set of modulated coefficients, before the step of checking of a detectability 
criterion. 

[0049] It is known to apply a predetermined distortion to a set of modulated coefficients and to carry out the check 
on the detectability criterion after this pre-distortion, in such a way that it is possible to guarantee the detection of the 
watermark inserted on a partition of the coefficients for a predetermined distortion. 

[0050] By virtue of the invention, a pre-distortion can be applied to the set of modulated coefficients, in a unique way, 
and does not need to be reiterated at each new region tested in the partition. 

[0051] In one particularly practical embodiment of the invention, the set of coefficients is a set of spectral coefficients 
originating from a spectral transformation of a digital image, and the step of applying a distortion includes the following 
sub-steps: 

inverse spectral transformation of the said set of modulated coefficients in order to obtain a watermarked image; 
application of a distortion to the said watermarked image; and 

spectral transformation of the watermarked image in order to obtain a set of modulated coefficients after distortion. 

[0052] The saving obtained in terms of complexity of calculation is all greater in this case whenever the application 
of a distortion requires several steps of processing of the set of coefficients, and particularly the application of an 
inverse spectral transformation in order to obtain the watermarked image and to apply a distortion, such as compres- 
sion/decompression to it, and then the application of a forward spectral transformation of the image again so as to 
check the criterion of detectability in the spectral domain. 

[0053] According to another particularly practical characteristic of the invention, the centred pseudo-random se- 
quence is formed from pseudo-random sub-sequences of identical size. 

[0054] Such a pseudo-random sequence is particularly well adapted to the application of a bottom-up or top-down 
iterative partitioning of the quadtree type which makes it possible iteratively to combine regions of identical sizes or 
iteratively to divide regions into sub-regions of identical sizes. 

[0055] For preference, the pseudo-random sub-sequences have a size greater than or equal to a minimum number 
of coefficients forming a region of the said partition of a size which is statistically significant for validating the said 
criterion of detectability over the said region. 

[0056] It is thus possible theoretically to form a maximum partition of the set of coefficients with the same unique 
centred pseudo-random sequence, each pseudo-random sub-sequence being centred on the regions with the smallest 
size which is statistically admissible for inserting a watermark by modulation into the set of coefficients. 
[0057] According to one practical embodiment, which is inexpensive in terms of calculations, at the generation step, 
each centred pseudo-random sub-sequence is created by generating half of the pseudo-random samples of the said 
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sub-sequence via a pseudo- random number generator and by symmetrisinq the said pseudo-random samples gener- 
ated 



[0058] For preference, in this embodiment, at the generation step, each centred pseudo-random sub-sequence is 
created by moreover permutating the said pseudo-random samples using a secret key. 

[0059] This permutation operation makes it possible better to preserve the intrinsic random-sequence nature of a 
sub-sequence created by symmetrising pseudo-random samples, and not to diminish the imperceptibility of the inserted 
watermark. 

[0060] According to another practical embodiment of the invention, at the generation step, each centred pseudo- 
random sub-sequence is created by generating pseudo-random random samples via a pseudo-random number gen- 
erator and by redistributing over each sample the sum of the said pseudo-random samples generated. 
[0061] According to a preferred characteristic of the invention, a partitioning by iterative division of the digital image 
is performed and, for each region at a rank of said partitioning, a partitioning of a directly higher rank of said region is 
validated if and only if there exist at least two subregions of said region for which said detectability criterion is validated. 
[0062] Thus, when it is attempted to divide a region into subregions, this partitioning of a higher rank is validated 
only if the number of regions in which a watermarking signal can actually be inserted whilst complying with the detecta- 
bility criterion is actually increased. 

[0063] Otherwise, it is preferable to keep the region before division, of greater size, in order to increase the detecta- 
bility of the watermarking signal. 

[0064] According to a preferred embodiment of the invention, for each region at a rank of said partitioning, the value 
of a capacity of the digital image is incremented when the partitioning of directly higher rank is not validated, the said 
capacity being equal to the cardinal number of the distinct regions obtained by the said partitioning on which the de- 
tectability criterion is validated. 

[0065] It is thus possible to calculate the capacity of the digital image which thus corresponds to the number of 
information bits which can be inserted in the digital image. 

[0066] By virtue of the adaptive partitioning performed by the invention, by reiterating the partitioning of each region 
as long as the detectability criterion is complied with, there is obtained the maximum capacity of the digital image which 
corresponds to the watermarking capacity of this image by the technique of insertion by modulation. This watermarking 
capacity can be defined as the maximum number of information bits which can be inserted imperceptibly in the imaqe 
with a guarantee of detectability on decoding. 

[0067] According to another embodiment of the invention, for each region at a rank of said partitioning, the value of 
a capacity is incremented when the partitioning of directly higher rank is not validated, said value of the capacity is 
compared with a fixed capacity value and the partitioning is reiterated for a region having a detectability amplitude 
which is the highest amongst all the other regions to be processed when said capacity value is less than said fixed 
capacity value. 

[0068] This partitioning method makes it possible, in this embodiment, to find the best possible partitioning in terms 
of detectability when a predetermined number of information bits have to be inserted, equal to said fixed capacity value. 
[0069] By selecting, during the partitioning, the regions for which the detection amplitude is the highest, the present 
invention makes it possible to determine the partitioning which, a priori, guarantees the best detectability on decoding 
of the watermarking signals making it possible to insert a predetermined number of information bits. 
[0070] According to another embodiment of the invention, a partitioning is effected by iterative merging of the digital 
image and, for each region of said partitioning for which said detectability criterion is not validated, said region is merged 
with another region of said partitioning. 

[0071] Preferably, in this embodiment, said other region of the partitioning is if possible a region for which said de- 
tectability criterion is not validated. 

« [0072] It is thus possible to increase the capacity of the digital image by joining the regions in which the detectability 
criterion is not satisfied. 

[0073] According to another aspect of the invention, a method of inserting, into a digital image, a watermark com- 
prising one or more information bits, comprises a method of determining a partition in accordance with the first aspect 
of the invention, and a step of inserting the information bit or bits by modulation of the coefficients of one or more 
50 regions of the said partition. 

[0074] This insertion method exhibits advantages similar to those of the method of determining a partition which it 
implements, and in particular makes it possible to insert a watermark in an adaptive and reliable way as reqards its 
detection. 

[0075] The present invention also concerns a method of decoding a watermarking signal in a digital image, inserted 
55 by modulation of coefficients representing said image in distinct regions forming a partitioning of said digital image. 
[0076] In accordance with the invention, this decoding method includes a step of determining a partitioning into 
distinct regions of the digital image to be decoded by effecting an adaptive partitioning as a function of a criterion of 
detectability of an information bit obtained by demodulation of the coefficients on each region. 
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[0077] Thus this decoding method makes it possible to find, at the decoder, the adaptive partitioning performed for 
inserting a watermarking signal. 

[0078] This is because, when the partitioning has been performed in an adaptive fashion as described above by a 
method of determining a partitioning in accordance with the first aspect of the invention, this partitioning, which depends 
on the content of the image, is not known to the decoder 

[0079] The decoding method according to the invention then consists of testing all the acceptable partitionings and 
verifying the detection of a watermarking signal on each of the regions associated with this partitioning. 
[0080] According to a preferred characteristic of the invention, the decoding method includes a prior step of comparing 
the size of each region of the partitioning with a minimum size corresponding to the minimum size of a statistically 
significant sample for the detection of a watermarking signal inserted in said region. 

[0081] According to a practical embodiment of the decoding method according to the invention, for each region of 
the partitioning, a detection amplitude is calculated from said detectability criterion and said detectability criterion is 
validated by comparing said detection amplitude with a threshold value predetermined for decoding. 
[0082] The decoding method thus makes it possible to find the adaptive partitioning effected when the watermarking 
signal is inserted. 

[0083] According to advantageous characteristic of this embodiment in the invention, the threshold value predeter- 
mined for decoding is less than said predetermined threshold value used during the method of determining a partition- 

[0084] This is because, at decoding, the detection test is performed on coefficients which have a priori been modu- 
lated and possibly made noisy. These coefficients are different from the original coefficients, and the result of the 
detectability test at decoding will be different from that performed in order to determine the partitioning. By choosing 
different thresholds on decoding and coding, it is possible to take into account the tendency for variation in detection 
amplitude due to an increased variance in the coefficients which have been modulated and possibly made noisy. 
[0085] Correlatively, the present invention also concerns a device for determining a partitioning into distinct regions 
of a digital image for inserting a watermarking signal by modulating coefficients of said regions, comprising partitioning 
means adapted to effect an adaptative partitioning as a function of a criterion of detectability of an information bit 
inserted on each region. 

[0086] According to a preferred embodiment, this device for determining a partition includes: 

• means for generating a centred pseudo-random sequence equal in size to the cardinal number of the said set of 
coefficients, formed by centred pseudo-random sub-sequences; 

- means for modulation of a set of coefficients by the said centred pseudo-random sequence in order to insert the 
same information bit on the said set of coefficients; and 

- means for checking, over each region, a criterion of detectability of the information bit inserted by modulation. 

[0087] The present invention also envisages a device for inserting into a digital image a watermark comprising one 
or more information bits, comprising a device for determining a partition in accordance with the invention, and means 
for inserting the information bit or bits by modulation of the coefficients of one or more regions of the said partition. 
[0088] Correlatively, the present invention also concerns a decoding device adapted to implement the decoding 
method according to the invention. 

[0089] The present invention also concerns a method of processing an image, comprising the following steps: 
inputting an image; 

inserting a first watermarking information into said inputted image; 
applying a distortion to the image obtained in said insertion step; 

- detecting said first watermarking information from the image obtained in said distortion applying step; and 

- determining an amount of a second watermarking information which can be inserted into said inputted image in 
accordance with the result of the detection obtained in said detection step. 

[0090] It further comprises the step of displaying the amount of the watermarking information determined in said 
determining step. 

[0091] Correlatively, the present invention also concerns an image processing apparatus, comprising: 
means for inputting an image; 

means for inserting a first watermarking information into said inputted image; 
means for applying a distortion to the image obtained by said insertion means; 

- means for detecting said first watermarking information from the image obtained by said distortion applying means; 
and 
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means tor determining an amount of a second watermarking information which can be inserted into said inputted 
image in accordance with the result of the detection obtained by said detection means. 

[0092] This apparatus further comprises the means for displaying the amount of the watermarking information de- 
termined by said determining means. 

[0093] The present invention also concerns a method of inserting a watermarking signal which can take into account 
at the same time predetermined imperceptibility and detectability constraints in order to ensure invisibility of the wa- 
termarking signal in the image and correct detection at the time of decoding. 

[0094] In the article entitled "A method for signature casting on digital images 0 by I. PITAS, in Proc. ICIP, pages 

215-218, September 1996, there is described a method of inserting a watermarking signal in a digital image in which 

the weighting coefficient a is determined so as to ensure detectability of the signal with a fixed probability. 

[0095] This weighting coefficient must be greater than a minimum value, which can be referred to as the detection 

amplitude, so as to allow the detection of the inserted signal with the correct detection probability level. 

[0096] This minimum value or detection amplitude depends notably on the length of the watermarking signal and on 

the required detection probability level. 

[0097] However, the weighting coefficient is fixed in the aforementioned document without taking into account the 
concept of psychovisual amplitude which guarantees invisibility of the watermarking signal inserted in the digital image. 
[0098] To this end the present invention relates to a method of inserting a watermarking signal in a set of coefficients 
representing a digital image, in which at least one subset of coefficients is modulated according to a linear model using 
a weighting coefficient, the weighting coefficient being greater than a minimum value determined as a function notably 
of the length of the watermarking signal and a detection probability level. 

[0099] According to the invention, this insertion method includes a step of calculating a maximum value of the weight- 
ing coefficient according to the length of the watermarking signal ensuring imperceptibility of the watermarking signal. 
[0100] Calculating a maximum value for the weighting coefficient, which can be referred to as visual amplitude, 
makes it possible to take into account, when a watermarking signal is inserted in a set of coefficients representing a 
digital image, constraints related to the imperceptibility of the inserted watermarking signal. 

[0101] Such an insertion method makes it possible to determine the conditions which make possible the insertion of 
a watermarking signal in an image whilst meeting predetermined imperceptibility and detectability criterion or on the 
other hand to decide if, under certain conditions, the insertion of the watermarking signal meets the predetermined 
imperceptibility and detectability criterion. 

[0102] According to a practical characteristic of implementation of the present invention, the insertion method also 
includes a step of determining compatible values of the length of the watermarking signal, of the weighting coefficient 
and of the probability level so that the value of the weighting coefficient determined is less than or equal to said maximum 
value and greater than or equal to said minimum value for said determined values of the length of the watermarking 
signal and of the detection probability level. 

[0103] By satisfying the aforementioned inequality for the weighting coefficient, the conditions of detectability and 
imperceptibility are satisfied at the same time. 

[01 04] It is thus possible to determine, in a correlated fashion, the value of the weighting coefficient, of the length of 
the watermarking signal and of the detection probability level in order to satisfy this inequality 

[01 05] According to a preferred characteristic of the invention, during the step of calculating the said maximum value, 
a function independent of said set of coefficients to be watermarked is used. 

[0106] Thus the model used for calculating the maximum value or the visual amplitude is independent of the image 
to be watermarked. 

[0107] In a preferred embodiment of the invention, the coefficients are spatio-frequency coefficients obtained by 
spatio-frequency transformation of this image, the calculation of said maximum value of the weighting coefficient de- 
pending on the type of transformation used. 

[0108] It is normal, before the insertion proper of a watermarking signal to apply a spatio-frequency transformation 
to the image, also referred to as transformation into sub-bands of the image, of the discrete cosine transformation 
(OCT) or wavelet transformation type, so that the modulation for inserting the watermarking signal is implemented on 
spatio-frequency coefficients of a sub-band of the image. 

[0109] The insertion method according to the invention applies particularly well to the conventional techniques of 
insertion by modulation inserted in a transformed domain obtained by one of the aforementioned spatio-frequency 
transformations. 

[0110] In practice, said calculation function used for calculating the maximum value or visual amplitude c^asa 
function of the length P of the watermarking signal w is of the type: 
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where 



a base( s ) is a base value, dependent on the transformation used S and on the (base) sub-band considered for the 
insertion, of the maximum weighting coefficient ensuring imperceptibility when a single coefficient of said sub-band 
is modulated, 

P is strictly greater than 2, and 

E[lwlP] is the mathematical expectation of the function IvvlP. 

[0111] According to another preferred characteristic of the invention, said minimum value is determined by using a 
calculation function also depending on the variance of the coefficients to be modulated, on the distribution of the wa- 
termarking signal and possibly on the variance of an additive noise decorrelated from the coefficients modelling any 
distortions applied to coefficients. 

[0112] Thus the minimum value of the weighting coefficient or detection amplitude can be calculated by taking into 
account notably the distortions which the image may undergo, notably during normal processings consisting of com- 
pression of digital images, in so far as such distortions can be modelled in the form of an additive noise decorrelated 
from the coefficients representing the image. 

[0113] In practice, the said function of calculating the detection amplitude a D as a function of the length P of the 
watermarking signal w is of the type 



a{Jo 2 x + c 2 n ) 



JbP^c 

where a* is equal to the variance of the coefficients to be modulated, a* is equal to the variance of the additive 
noise and a, b and c are constants which depend on the distribution of the watermarking signal and on the required 
detection probability level. 

[0114] According to a preferred characteristic of the invention, the watermarking signal is a pseudo-random sequence 
with a predetermined distribution with a zero mean. 

[0115] The use of such a watermarking signal increases the robustness of the inserted watermarking signal with 
respect to an intentional attack aimed at extracting or modifying the watermarking signal, typically by statistical analysis, 
the kernel of the pseudo-random sequence produced from a predetermined distribution also being able to be qiven by 
a secret key. 

[0116] In a particularly advantageous embodiment of the invention, the method of inserting a watermarking siqnal 
includes the following steps: 

- calculating a maximum value of the weighting coefficients for a length of the watermarking signal equal to the 
cardinal number of the set of modulatable coefficients ensuring imperceptibility of the watermarking signal at the 
limit of perceptibility; 

- calculating the length of the watermarking signal for a minimal value of the weighting coefficient equal to said 
maximum value calculated and for a predetermined detection probability level; and 

- comparing said calculated length with the cardinal number of the set of modulatable coefficients. 

[0117] Thus, by determining the maximum value of the weighting coefficient at the limit of perceptibility, that is to say 
in the eventuality of all the coefficients available being modulated, it is possible to determine the minimum length of 
the watermarking signal which makes it possible to obtain a correct detection probability with a weighting coefficient 
equal to the maximum value calculated. 

[01 1 8] The comparison of this length thus calculated with the cardinal number of all coefficients available in order to 
effect the insertion makes it possible to determine whether this insertion is actually possible with the required detection 
probability level. 

[0119] In addition, if the calculated length is less than the cardinal number of the set of coefficients, it can be certain 
that the imperceptibility criterion will indeed be met by modulating a subset of coefficients with a cardinal number equal 
to this calculated length, since the value of the weighting coefficient was calculated in the eventuality of all the available 
coefficients actually being modulated. 

[0120] In practical terms, if, at said comparison step, said calculated length is less than or equal to the cardinal 
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number of the set of modulatable coefficients, a watermarking signal with a length al ieast equal to said calculated 
length is inserted by modulating a subset of coefficients with a cardinal number equal to said length, according to a 
linear model using a weighting coefficient equal to said calculated maximum value. 

[0121] Conversely, if at said comparison step said calculated length is greater than the cardinal number of the set 
of modulatable coefficients, there is calculated the detection probability level obtained for a minimum value of the 
weighting coefficient equal to said maximum value calculated for a watermarking signal length equal to the cardinal 
number of all the modulatable coefficients. 

[0122] Thus the decision to insert a watermarking signal or not on the set of modulatable coefficients can be taken 
in the theoretical knowledge of the detection probability level obtained by default. 

[0123] In another practical embodiment of the invention, the insertion method includes a step of calculating a thresh- 
old value of the length of the watermarking signal determined so that said minimum value of the weighting coefficients 
is equal to said maximum value of the weighting coefficients for a length of the watermarking signal equal to said 
threshold value and a predetermined detection probability level. 

[0124] It is thus possible to determine the minimum length threshold which the watermarking signal must meet in 
order to simultaneously meet the two criteria of perceptibility and detection when the watermarking signal is inserted 
in a digital image. 

[0125] Correlative^, the present invention also concerns a device for inserting a watermarking signal in a set of 
coefficients representing a digital image, in which at least one subset of coefficients is modulated according to a linear 
model using a weighting coefficient, the weighting coefficient being greater than a minimum value determined as a 
function notably of the length of the watermarking signal and of a detection probability level. 

[0126] This insertion device comprises means of calculating a maximum value of the weighting coefficient as a func- 
tion of the length of the watermarking signal ensuring imperceptibility of the watermarking signal. 
[0127] This insertion device has characteristics and advantages similar to those described previously since it is 
adapted to implement the insertion method according to the invention. 

[0128] The present invention also concerns a computer, a digital image processing apparatus, a digital printer, a 
digital photographic apparatus and a digital camera adapted to implement the method of determining a partitioning, 
the methods of inserting a watermark, the method of processing an image and the decoding method in accordance 
with the invention. 

[0129] These devices for determining a partitioning, inserting a watermark, processing an image and decoding, this 
computer, this digital image processing apparatus, this digital printer, this digital photographic apparatus and this digital 
camera have characteristics and advantages similar to those described with reference to the methods which they 
implement. 

[01 30] The present invention also relates to a computer program product and to an information storage means which 
can be read by a computer or a microprocessor, possibly removable, which stores a program implementing the method 
of determining a partitioning and/or inserting a watermark and/or processing an image and/or decoding in accordance 
with the invention. 

[0131] Other particularities and advantages of the invention will also emerge from the following description of the 
invention. 

[0132] In the accompanying drawings, given by way of non-limitative examples: 

Figure 1 is an overall diagram illustrating a device for inserting a watermark into a digital image; 
Figure 2 schematically illustrates a step of spectral decomposition employed when inserting a watermark- 
Figure 3 is an algorithm of the method of determining a partitioning in accordance with a first embodiment of the 
invention; 

Figure 4 is an example of partitioning of a region by division; 

Figure 5 is an algorithm of the method of determining a partitioning in accordance with a second embodiment of 
the invention; 

Figure 6 is an algorithm of the decoding method according to a first embodiment of the invention; 

Figure 7 is an algorithm illustrating the application of the method of determining a partitioning of Figure 1 to a digital 

image; 

Figure 8 is an algorithm illustrating the application of the decoding method of Figure 4 to a digital image; 
Figure 9 is a block diagram illustrating a variant embodiment of the method of determining a partitioning in accord- 
ance with a third embodiment of the invention; 
- Figure 10 illustrates the application of the method of determining a partition to a digital image according to a 
preferred embodiment of the invention; 

Figure 11 is an algorithm of the method of determining a partition in accordance with a first way of the preferred 
embodiment of the invention; 

Figure 1 2 is an algorithm of the method of determining a partition in accordance with a second way of the preferred 
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embodiment of the invention ; 

- Figure 13 is an algorithm of the insertion method according to one embodiment of the invention; 

- Figure 14 illustrates, by means cf a graph, the insertion method of the figure 13; and 

• Figure 1 5 is a block diagram illustrating a device adapted to implement the methods in accordance with the inven- 
tion. 

[0133] A description will first of all be given, by reference to Figure 1, of the insertion of a watermark into a set oi 
coefficients representative of a digital image I. 

[01 34] This digital image I can be represented by a set of coefficients, either in the spatial domain or in a transformed 
domain, the coefficients in this latter case being hybrid coefficients, that is to say coefficients located both in the spatial 
domain and in the frequency domain. Such a representation of the image is obtained, for example, by using a sub- 
band originating from a space-frequency decomposition of the image, for example a discrete wavelet decomposition 
[0135] Here, the insertion of the watermark is achieved via a spread-spectrum insertion technique by modulation of 
coefficients of a space-frequency representation of the image, obtained by the space-frequency transformation of the 
digital image I. 

[01 36] By way of example, a discrete wavelet decomposition S can be used, represented diagrammatically in Figure 
2. This space-frequency decomposition is well known in the image-processing field, and only the principle thereof is 
reiterated below. It makes it possible to segment the image into sub-bands of frequencies and to obtain hybrid coeffi- 
cients, that is to say spectral coefficients which are also located in space, here in the plane of the image 
[0137] The image I consists of a series of digital samples. The image I is represented, for example, by a series of 
bytes, each byte value representing one pixel of the image I , which may be a black-and-white image with 256 grey levels 
[0138] The multi-resolution spectral decomposition means consist of a sub-band decomposition circuit or analysis 
circuit, formed by a set of analysis filters, associated respectively with dividers-by-two. This decomposition circuit filters 
the image signal I in two directions, into sub bands of high spatial frequencies and of low frequencies The circuit 
includes several successive analysis blocks for breaking down the image I into sub-bands according to several levels 
of resolution. 

[01 39] By way of example, the image I is broken down here into sub-bands with a level of decomposition d equal to 3 
[01 40] A first analysis block receives the image signal I and filters it through two respectively low-pass and high-pass 
digital filters, in a first direction, for example horizontal. After passing through dividers-by-two, the resultant filtered 
signals are in their turn filtered by two respectively low-pass and high-pass filters, in a second direction, for example 
vertical. Each signal is again passed through a divider-by-two. Hence four sub-bands LL 1f LH 1( HL, and HH are 
obtained at the output of this first analysis block, with the highest resolution in the decomposition.' 1 
[0141] The sub-band LL, includes the low-frequency components in the two directions of the image signal I The 
sub-band LH 1 includes the low-frequency components along a first direction and high-frequency components along a 
second direction of the image signal I. The sub-band HL, includes the high-frequency components along the first 
direction and the low-frequency components along the second direction. Finally, the sub-band HK includes the high- 
frequency components along the two directions. 

[0142] A second analysis block in its turn filters the sub-band of low frequencies LL, so as, in the same way, to supply 
four sub-bands LL 2 , LH 2 , HL 2 and HH 2 with intermediate level of resolution in the decomposition. A third analysis block 
then filters the sub-band of low frequencies LL 2 so as to supply four sub-bands LL 3 , LH 3 , HL 3 and HH 3 . Finally in this 
example, the sub-band LL 3 is in its turn analysed by a fourth analysis unit in order to supply four sub-bands LL 4 , LH 4 , 
HL 4 and HH 4 with the lowest resolution in this decomposition. 

[0143] Hence thirteen sub-bands and four levels of resolution are obtained. Obviously, the number of levels of res- 
olution, and consequently of sub-bands, may be chosen differently. 

[0144] Obviously, other types of spectral transformation could be used, such as the discrete Fourier transform the 
discrete cosine transform or the Fourier-Mellin transform. 

[0145] In a general way, frequency sub-bands are obtained forming a set of spectral coefficients into which a water- 
mark can be inserted. 

[01 1 46] The high-frequency sub-band of highest resolution HH, can be considered here in order to insert a watermark. 
This sub-band HH, thus supplies a set of coefficients X, equal to N in size, for example. 
[0147] Let X = fx,, 1 < i < N) be the set of coefficients representative of the digital image I. 

[0148] Obviously, it would be possible to consider a subset of coefficients of this set X, representative only of a part 
of the digital image I. 

[0149] In order to insert a watermark into this set of coefficients X, one technique consists in inserting a pseudo- 
random signal, by spreading its spectrum, so as to render this signal undetectable by spectral or statistical analysis., 
[0150] Let w = {Wj, 1 < j < P} be a watermark of size P < N, a pseudo-random signal with known distribution and 
average of zero. The most widespread distributions are the binary distribution {-1 , 1 }, the uniform distribution over [-1 
1] and a centred normalised Gaussian distribution N (0, 1 ). 
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[0151] The linear insertion formula is : 

XV = Xj + b a. w ]( with 1 < j < P, 

in which {Xj, 1 < j < P} is a subset of the set of coefficients X, b is an information bit, and ctj is a weighting coefficient, 
also called modulation amplitude. 

[0152] The set {j, 1 < j < P} is also called information-bit insertion carrier. 

[0153] In order to insert a signal which may be formed from several information bits, it is worth partitioning the set 
of coefficients so as to determine the number of regions or the number of insertion carriers available in this set. 
[0154] To do that, the insertion device 10 in accordance with the invention comprises a device for determining a 
partition 11 and insertion means 12. 

[0155] The device for determining a partition 11 is suitable for carrying out adaptive partitioning of the set of coeffi- 
cients as a function of the detectability of an inserted watermark. 

[0156] This method of adaptive partitioning consists in varying the number of coefficients modulated in order to insert 
an information bit of the watermark, so as simultaneously to satisfy the criteria which are set of invisibility and of correct 
detection probability. 

[0157] A description will first of all be given, with reference to Figure 3, of a method of determining a partitioning into 
distinct regions of a digital image which makes it possible to maximise the number of regions obtained in order to 
maximise the number of information bits inserted in the image. 

[0158] Consideration is given to a representation of the digital image by a set of coefficients, either in the spatial 
domain, or in a transformed domain, the coefficients being in the latter case hybrid, that is to say located both in the 
spatial domain and in the frequency domain. Such a representation of the image is for example obtained by using a 
sub-band issuing from a spatio-frequency decomposition of the image, as a discrete wavelet decomposition previously 
described. 

[0159] The partitioning of the image is effected in order to insert a watermarking signal by modulation of the coeffi- 
cients of each region. 

[01 60] According to the invention, this partitioning into regions is performed by an adaptive partitioning as a function 
of a criterion of detectability of an information bit inserted on this region. 

[0161] The number of regions obtained under constraint of the detectability criterion makes it possible to define a 
maximum capacity of the digital image. This maximum capacity corresponds to the number of information bits which 
it is possible to insert in the digital image whilst complying with the detectability criterion. 

[0162] Initially, at a step E100, this capacity Q is initialised to 0 and the starting point is an initial partitioning which 
is limited here to a single region R of size P. 

[0163] Naturally, the initial partitioning could be different and notably include already several distinct regions with 
sizes which might be different. 

[0164] It is assumed that the size P of this initial region R is greater than a minimum size denoted P slat which corre- 
sponds to the minimum size of a statistically significant sample of coefficients for detecting a watermarking signal 
inserted in these coefficients. Typically, this minimum size can be fixed at 100. . 

[0165] This initial region R can consist of the set of coefficients of the representation under consideration or a subset 
of this representation. 

[0166] A test step E101 then tests the detectability of the insertion of a watermarking signal in the initial region R. 
[0167] This detectability test can be effected by actual modulation and detection operations on the coefficients of 
the initial region R. 

[0168] In practice, all the coefficients of the initial region R are copied into a working memory so as not to directly 
watermark the image. 

[0169] Next all the coefficients of the initial region R are modulated by inserting the watermarking signal. 
[0170] Conventionally, the following procedure is followed in order to effect this modulation: 
[0171] The initial region R consists of a set of coefficients: 

X: {X |( i= 1, ...P} 
[0172] Normally a linear modulation is used of the type: 

X j * = X i ± (XjW., withi = 1 P 

where W : {Wj, i = 1 P} is the watermarking signal generally chosen as a pseudo-random sequence with a 

known distribution and a nil mean. 

[0173] The value of the sign ± of the modulation depends on the binary value to be inserted: for example the sign - 



11 



EP 1 043 687 A1 



corresponds to the value 0 and the sign + corresponds to the value 1 . 

[0174] It is possible to choose, by way of example, a pseudo-random sequence VV which follows a uniform law on 
the interval [-1, 1]. 

[0175] Naturally, any pseudo-random watermarking with known distribution and nil mean can be suitable. The most 
usual distributions for the watermarking signal W are, apart from the aforementioned uniform distribution [-1, 1), the 
binary distributionM ,1 } and the centred standardised Gaussian distribution N(0, 1). 

[0176] The modulation can optionally be protected by means of a secret key characterising for example the kernel 
for reproducing the pseudo-random sequence W. 

[0177] The term ctj is the modulation amplitude, or weighting coefficient, applied to the coefficient Xj. 
[0178] A constant weighting coefficient can be used for all the coefficients of the region R, such that c- = <x v for any 
i, the value of the weighting coefficients being equal to a maximum value ensuring imperceptibility of the watermark- 
ing signal for the modulation of the set of coefficients representing said digital image. 

[0179] This maximum value oc v (referred to as Just Noticeable Difference or JND amplitude) decreases in fact with 
the length of the sequence of the watermarking signal W, and consequently with the number of modu latable coefficients. 
[0180] Given that the capacity Q will be maximised when all the coefficients of the representation of the digital image 
are used for the modulation, it is judicious to determine a modulation JND amplitude for a number of modulatable 
coefficients equal to the cardinal number of the set of coefficients representing the image. 

[0181] This maximum value a v or JND amplitude corresponds to the maximum modulation amplitude which can be 
used in the aforementioned linear insertion model beyond which an observer is capable of visually detecting a change 
in the reconstructed watermarked image. 

[0182] In order to determine this JND amplitude, it is possible to use a visibility model which makes it possible to 
predict the visibility of a watermarking operation according to different parameters, namely the representation of the 
signal through a spatio-frequency transformation S used, the sub-band considered for the insertion, the type of distri- 
bution of the sequence of the watermarking signal W and the length P of the sequence W. 

[0183] A simple model, developed by WATSON and described in the article entitled "Visibility of wavelet quantization 
noise", A.B. WATSON et al, IEEE Trans, on Image Process, 6(B), 1164-1175, 1997, makes it possible to predict the 
visibility of a set of modulated coefficients from the measurement of visibility of a single modulated coefficient. Reference 
can advantageously be made to this document for the detailed description of this model. 

[0184] Thus use is made of a function dependent on the length P of the watermarking signal W, the type of trans- 
formation S used and the (base) sub-band considered for the insertion, but independent of all the coefficients to be 
watermarked X. 

[0185] This calculation function used can be written: 



a IP S w\ - a t>Bs* {S) 



where 



«base( s ) is a °ase value, dependent on the transformation used S and on the (base) sub-band considered for the 
insertion, of the maximum weighting coefficient ensuring imperceptibility during modulation of a single coefficient 
of this sub-band, 
P is strictly greater than 2, and 

E[\w\V] is the mathematical expectation of the function livll 3 . 



[0186] The base values a^^S) can be measured once and for all, for each sub-band of coefficients in the wavelet 

decomposition, from a single psychovisual measurement, and be stored in a table of visual amplitudes. 

[0187] p is the exponent of a Minkowsky sum and can be chosen equal to 5, for example. 

[0188] The mathematical expectation E[lw|P] corresponds to an estimation of the mean of the function MP 

[01 89] This visibility model does not take into account the image I to be watermarked itself and is independent thereof. 

This is equivalent to considering that the image I is uniform. It is a question of a "worst case" model since the presence 

of the image signal makes it possible to visually mask the modulation itself. 

[0190] The JND amplitude can thus be calculated for a length P of the watermarking signal W equal to the cardinal 
number of all the modulatable coefficients of the sub-band under consideration. 

[0191] However, it is judicious to exploit the fact that the image signal itself makes it possible to mask the modulation. 
[0192] Thus, in this example, use is made, for each coefficient Xj to be modulated, of a weighting coefficient in 
accordance with a law of the form = kj.cv, where k t is a modulation factor dependent on the coefficients situated 
close to the coefficient considered on the region R and a v is equal to the maximum JND value of the weighting coef- 
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[0193] Thus each coefficient will be modulated according to the local content, which makes it possible to locally 
increase the amplitude of the modulation to the benefit of detectability. 

[0194] In order to test detectability in the test step E101, a delegability amplitude is next calculated from the de- 
legability criterion and this detectability criterion is validated by comparing the detectability amplitude with a predeter- 
mined threshold value denoted Tc. 

[0195] A calculation of correlation C(X*. W) between all the modulated coefficients X* and the marking signal W is 
generally made and it is decided whether there has actually been watermarking of the image when the result of the 
correlation calculation C(X*. W) is greater than the predetermined threshold value Tc. 

[0196] A standardised statistical test can also be used for the detection, such as the one described in the article 
entitled "A method for signature casting on digital images" by I. PITAS, in Proc. ICIR pages 21 5-218, September 1 996 
The detection is then characterised in terms of probability. It is thus possible to choose a threshold value Tc corre- 
sponding to a fixed detection probability level, for example 99.95%. 

[0197] Thus, for a given region R defining a set of coefficients X of size P and with a fixed weighting coefficient a it 
is possible to calculate a detectability amplitude and compare it with a threshold value Tc so as to validate or not the 
detectability criterion on this region R. 

[0198] The detectability criterion is in fact in theory a function depending on the size P of the region R, the detection 
threshold Tc related to the required detection probability, the watermarking signal W, the variance o\ of the coefficient 
X to be modulated on the region R and the weighting coefficient used a for modulating the coefficients X according to 
the linear model described above. 

[0199] Thus, for example, for a detection by standardised hypothesis test and a value of the weighting coefficient a 
which is constant and equal to a v for all the coefficients Xj of the region R, a minimum length P min (X) is defined by the 
function: 

Pmin(X)= (^) 2 (^ 2 x+ ca + d\ 

where a, b, c and d can be assimilated to constants. The signal being non-stationary, that is to say the variance 
o2 x being a function of space, it is indeed necessary to adapt the modulation length P generated by the watermarkinq 
signal to the image signal itself. 

[0200] On decoding, the set of coefficients X being modulated, this function also depends on the variance a 2 w of 
the inserted watermarking signal W. 

[0201] In addition, it is necessary to provide for any distortions which the image may undergo after modulation, due 
for example to a compression of the image. 

[0202] When this noise can be considered to be additive and decollated from the coefficients, it is possible to use 
a model dependent on the variance o\ of this additive noise modelling any distortions applied to said coefficients after 
modulation. 

[0203] Thus, to a given weighting coefficient a, making it possible to ensure invisibility of the inserted watermarking 
signal, a detection threshold Td fixed as a function of the detection probability required at decoding and a set of selected 
coefficients X, there also corresponds a minimum modulation length P mjn (X) which makes it possible to ensure de- 
tectability of the inserted watermarking signal. 

[0204] On decoding, for a detection by standardised hypothesis test, this minimum length P min (X) is then defined by 
the function: y 

P min(X) = (^cT T [b{ ° 2 * + a2 °w + °*> + c*+ct) 

where a, b, c and d can be assimilated to constants. 
[0205] It is therefore necessary, if the same detectability test is used during coding and decoding, to choose a thresh- 
old value Tc on coding greater than that Td used at the time of decoding if it is wished to be able to find the same 
partitioning on decoding of the image. 

[0206] In practice, the threshold value Td on decoding which it is wished to use is fixed, that is to say for example 
the required detection probability level when the detection is made by a standardised hypothesis test, and a higher 
threshold Tc on coding is chosen, in order to take account of the variance in the watermarking signal W and the variance 
of an additive noise fixed a priori so as to anticipate the influence of the subsequent processing undergone by the image. 
[0207] This condition thus makes it possible to find the same detectability criterion which constrains the partitioning 
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of the digital image without having to transmit to the decoder the partitioning effected at the time of insertion of the 
watermarking signal W. 

[0208] In this example embodiment, a partitioning by iterative division of the digital image is effected. It is conse- 
quently considered, at step E101 , that the detectability criterion is satisfied when the calculated detection amplitude is 
greater than or equal to the predetermined threshold Tc. 

[0209] If the detectability criterion is not satisfied, then it is probable that the detectabiiiiy criterion is also not satisfied 
in subregions of the initial region R of smaller size because the detection probability generally decreases with the size 
of the regions. A step E 102 of stopping the method of determining a partitioning makes it possible to stop the process 
of iterative division of the image. 

[0210] Insertion not being judged possible in the initial region R, the capacity Q of the image is therefore nil. 
[0211] If on the other hand the detectability criterion is satisfied, that is to say here that the detection amplitude is 
greater than or equal to the threshold value Tc, it is possible that it is also satisfied for subregions of the initial region R. 
[0212] It is therefore judicious to attempt to divide the region R as long as the detectability criterion is satisfied in 
order to increase the capacity Q. 

[0213] There are of course several possible means for effecting the spatial division of the region R into subregions 
R j: of sizes Pj, with j = 1 K. 

[0214] In this example a division into a quadtree is chosen, as illustrated in Figure 4. 
[0215] This division consists of dividing each region into four subregions of identical size. 

[0216] Thus, by way of example, in Figure 2, a digital image of initial size N is being partitioned into two sub-regions 
R5 and R6 of sizes N/4, seven subregions R1, R2, R3, R4, R7, R8, R9 of sizes N/16 and four subregions R10 R11 
R 1 2, R1 3 of sizes N/64. 

[0217] Thus, at the division step E103, the region R is divided into subregions R, j = 1 K of sizes P.. In a case of 

a division into a quadtree, K = 4 and Pj = P/4. 

[0218] At step E104 a first subregion of index j (for example j = 1) is considered and a value of potential gain in 
capacity denoted G Q is initialised to - 1 . 

[0219] The partitioning determination method then includes a prior step E105 of comparing the size P } of the region 
Rj with the minimum size P stat corresponding to the minimum size of a statistically significant sample for the detection 
of a watermarking signal inserted in said region as described above. 

[0220] If the size Pj of the region R } is less than this minimum size P stat , at step E1 08 the index J = J+1 is incremented 
in order to process the following sub-region R^,. 

[0221] Otherwise, the detectability criterion is checked on the subregion R in the same way as previously described 
with reference to step E101 . 

[0222] In practice, the watermarking signal is inserted on the coefficients of the said region R jf a detection amplitude 
is calculated and this detectability amplitude is compared with the threshold value Tc. 

[0223] The detection amplitude can be a correlation value between the modulation signal X* and the carrier of the 
watermarking signal W or the result of a standardised hypothesis test. 

[0224] If the detectability criterion is verified, the value of potential gain in capacity G Q is increased by one unit at 
step E 107. 

[0225] Then at step E108 the index j of the subregions is incremented and it is checked at the test step E109 that 
this index j is less than or equal to the number of K. 

[0226] In the affirmative, steps E105 to E109 are reiterated on the following subregion. 

[0227] Otherwise, when the index j = K+1, it is checked in a test step E110, whether the potential gain in capacity 
value G Q is strictly positive, that is to say whether there exist at least two subregions R of the initial region R for which 
the detectability criterion is validated. 

[0228] In the negative, the partitioning of higher rank carried out on the initial region R is not validated and this region 
R is preserved. Such is the case when the detectability criterion has not been validated for any of the subregions R. 
or only for a single subregion Rj. 1 
[0229] In the latter case, it is preferable to keep the region R for inserting the watermarking signal since it is of size 
P greater than the size of the single validated subregion Rj. 

[0230] At step E112, the value of the capacity Q is then increased by one unit and the region R of size P in the 
partitioning is kept. 

[0231] An actual insertion step E1 1 3 by modulation of the coefficients of the region R can then be implemented in a 
normal fashion using the linear insertion model described above. 

[0232] An elimination step El 14 then makes it possible to eliminate the region R from a stack of regions to be proc- 
essed. 



[0233] On the other hand, if at step E1 10 the value of the potential gain capacity Gc is strictly positive, in an addition 
step E 1 1 1 all the subregions Rj for which the detectability criterion is validated in the stack of regions R to be processed 
are added. 
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[0234] A return step E115 makes it possible to reiterate all the steps E103 to E113 in order to process the regions 
of the stack in a recursive fashion. 

[0235] At the end of the recursive partitioning method just described, a capacity value Qm corresponding to the 
maximum capacity o! the image for insertion of a watermarking signal is obtained, guaranteeing a correct detection of 
this signal. This maximum capacity Qm is equal to the cardinal number of the distinct regions obtained by the said 
partitioning on which the detectability criterion is validated. 

[0236] In addition, the insertion step proper E113 can be omitted during a first "pass" which then makes it possible 
only to determine the maximum capacity Qm of the image and the validated regions. 

[0237] Once the maximum capacity Qm has been calculated and the corresponding partitioning is stored in memory, 
it is then possible for a user to determine the message which he wishes to insert, corresponding to a number of infor- 
mation bits equal to the maximum capacity Qm, and then to insert this message in a normal manner by modulating 
the coefficients of the different regions obtained by the previously determined partition. 

[0238] A description will now be given, with reference to Figure 5, of a second embodiment of the invention in which 
it is wished to insert a message of fixed capacity Qs in an image I. Naturally, this fixed capacity Qs must have a value 
less than the maximum capacity Qm of the digital image to be watermarked. 

[0239] As illustrated in Figure 5, the method of determining a partitioning in this embodiment includes steps E 100 to 
E1 13 identical to those of the method of determining a partitioning of the embodiment of Figure 3 and which it is not 
necessary to describe again here. 

[0240] As before, at step E 1 1 2, the value of the capacity Q is incremented when the partitioning of rank directly higher 
is not validated. 

[0241] Then the method includes a comparison step E116 in which this value of the capacity Q is compared with the 
fixed capacity value Qs. 

[0242] In practice, it is tested at step E116 whether the capacity value Q is equal to the fixed capacity value Qs. 
[0243] In the affirmative, a stop step E1 17 makes it possible to stop the iterative partitioning of the image since there 
are a number of regions equal to the number of information bits Qs which it is wished to insert. 

[0244] Otherwise, in an elimination step E1 14, the processed region is eliminated from the stack of regions R to be 
processed formed at step El 11. 

[0245] The method of determining a partitioning also includes, in this embodiment, a sorting step E118 which sorts 
the regions R to be processed in the stack by detectability amplitude value calculated on these regions R. 
[0246] Thus a return step E115 reiterates the partitioning from the division step E103 for a region R in the stack 
having the highest detectability amplitude amongst all the other regions to be processed. 

[0247] A sorting variant can consist of sorting the regions to be processed in decreasing order of size so as to process, 
as a priority, the largest regions which a priori provide better detectability. 

[0248] Naturally, a partitioning by iterative merging of the digital image could be carried out in the same way, for 
example also using a quadtree partitioning structure. 

[0249] The initial partitioning then consists of regions with smaller acceptable sizes, greater than or equal to the 
previously defined minimum size P stat . 

[0250] For each region of the partitioning for which the delectability criterion is not validated, this region is merged 
with another region of the partitioning. 

[0251] This other region of the partitioning is if possible a region for which the detectability criterion is not validated 
so as to increase the capacity 

[0252] Alternatively, another region can be a region satisfying the detectability criterion so as to make detection 
reliable. 

[0253] There are many possibilities for the choice of the region to be merged. It is possible to take as the merger 
criterion the detection amplitude calculated on each of the regions and merge for example the regions associated with 
the lowest detection amplitudes. 

[0254] In this particular case of a merging in a quadtree partitioning structure, the merging is structurally constrained 
and consists of deciding whether or not four adjacent blocks, for example R1, R2, R3 and R4 in Figure 4, must be 
merged. 

[0255] A description will now be given, with reference to Figure 6, of the associated decoding method which makes 
it possible to find the partitioning when the watermarking signal inserted in a digital image is extracted. 
[0256] Indeed, because of the adaptation of the partitioning to the digital image, this partitioning will not be known a 
priori at the decoder. 

[0257] The decoding method then consists, in general terms, of testing all the acceptable partition in gs. The same 
iterative implementation technique is used as the one used when the partitioning is determined for inserting the wa- 
termarking signal. 

[0258] Thus a description will be given of an example embodiment in the case of a segmentation by dividing by 
blocks in a quadtree, symmetrical to that using the partitioning determination method of the embodiment described 
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with reference to Figure 3. 

[0259] At an initialisation step E200, there is an initial partitioning limited to a single region R of size R This region 
consists of all the coefficients of the same representation of the digital image as the one used at the time of insertion 
of the watermarking signal. 

[0260] A test step E201 makes it possible to verify that the size P of the initial region R is indeed greater than the 
minimum size P stat as described previously. 

[0261] In the affirmative, at step E203 a criterion of detectabiiity of the information bit inserted on the region R, 
obtained by demodulation of the coefficients on this region R, is verified. 

[0262] In practice, a detection amplitude is calculated for all the coefficients of the region R and the detectabiiity 
criterion is validated by comparing this detection amplitude with a predetermined threshold value Td for decoding. 
[0263] It is possible to proceed by means of a standardised statistical test, and the result compared with the prede- 
termined threshold value Td. 

[0264] As explained before, it is necessary, if the same detectabiiity test is used during coding and decoding, to 
choose a threshold value Tc on coding greater than that Td used at decoding if it is wished to be able to find the same 
partitioning on decoding of the image even when the latter has been made noisy by post-processing. 
[0265] This condition thus makes it possible to find the same detectabiiity criterion which constrains the partitioning 
of the digital image without having to transmit, to the decoder, the partitioning effected when the watermarking signal 
W was inserted. 

[0266] If the detectabiiity criterion is not verified at step E203, a division is carried out, at step E204, of the region R 

into subregions Rj and sizes Pj with j = 1, K. Then the same partitioning structure is used as the one used on coding, 

that is to say, in this example, a division of the region R into four subregions Rj of identical sizes. 

[0267] Then, in an addition step E205, the subregions Rj are added in a stack of regions R to be processed and, at 

a return step E202, the decoding method is reiterated recursively as from the test step E201 on each region R of the 

stack. 

[0268] When at the step of calculating the detectabiiity criterion E203, the latter is verified for a region R, there is 
extracted at the extraction step E206 the information bit inserted in this region R as an element of the method inserted 
in the digital image I. This extraction is effected in a conventional fashion by calculating a detection amplitude. The 
result of the test gives, in absolute value, a detection amplitude to be compared with the threshold value Td of the 
decoding, the sign giving the value 0 or 1 of the information bit inserted. 

[0269] It should be noted that this extraction is already partly effected during the detection step E203 so that the 
detection E203 and extraction E206 steps are structurally coupled. 

[0270] An elimination step E207 is implemented in order to eliminate the region R from the stack of regions to be 
processed, which was formed at the addition step E205, and the return step E202 then makes it possible to implement 
the decoding method recursively on each region R of the stack to be processed. 

[0271] A description will now be given, with reference to Figures 7 and 8, of a practical application of the method of 
determining a partitioning of a digital image in a particular embodiment of the invention. 

[0272] In this application, the watermarking signal is inserted by means of a technique of insertion by spectrum 
spreading, by modulation of coefficients of a space-frequency representation of the image, obtained by a spatio-fre- 
quency transformation S of the digital image I as previously described with reference to figure 2. 

[0273] In general terms Sj sub-bands are obtained, with J = 1 M and M equal here to 1 3 at a transformation step 

E300. 

[0274] Next, for each sub-band Sj of index j considered at step E301, a partitioning is effected at step E302 by 
implementing the method of determining a partitioning described previously to with reference to Figure 3. 
[0275] In this way a capacity Qj is derived for each sub-band Sj. 

[0276] A step E303 of insertion by modulation of the coefficients next makes it possible to insert Q: information bits 
in the Qj regions validated. 

[0277] In steps E304 and 305, the following sub-bands are next considered as long as all the sub-bands Sj have not 
been processed, and the steps E302 and E303 of partitioning and inserting a watermarking signal of capacity Q are 
reiterated on each sub-band Sj. 

[0278] Next, an inverse transformation S 1 is applied to the image in order to obtain the watermarked image, of total 
capacity Q equal to the sum of the capacities Qj calculated for each of the sub-bands Sj. 

[0279] Thus, each sub-band of the representation is considered independently for the adaptive partitioning which is 
carried out both at coding and at decoding. 

[0280] As illustrated in Figure 8, the decoding method in this embodiment is implemented in a similar fashion using 
the decoding method described with reference to Figure 6, the initial region R here corresponding each time to a sub- 
band of the decomposition of the image. 

[0281] A step E400, identical to step E300, divides the image I into Sj subbands with j = 1, ... M. 

[0282] A first sub-band Sj is considered at step E401 . The order of processing of the sub-bands is identical to that 



16 



EP 1 043 687 A1 



used during coding in order to extract the inserted message in order. 

[0283] The adaptive partitioning is implemented at the partitioning step E402 in order to calculate the capacity Q of 
the sub-band Sj. Next, at the extraction step E403, the Qj information bits are extracted by demodulating the coefficients 
of the Qj regions of the partitioning. 

[0284] At the test steps E404 and E405, it is checked in the usual fashion that all the sub-bands have been processed 
and, in the negative, the decoding steps E402 and E403 (partitioning and extraction) are reiterated on the remaining 
sub- bands. 

[0285] In this way, at the output of the decoder, a message of Q bits is obtained, Q being equal to the sum of the 
capacities Qj of each sub-band Sj. 

[0286] In a variant of this method illustrated in Figure 9, it is possible to effect a partitioning on a post-processed 
image signal in order to obtain an a priori robustness of the watermarking signal to any post-processing undergone by 
the image. 

[0287] The principle of this variant consists of performing a detectabiiity test in a test step E309 which is no longer 
coupled to the steps of partitioning E302 and insertion E303 described previously. 

[0288] After the insertion of the watermarking signal in a region, a reverse transformation S" 1 is first of all performed 
at step E306 in order to find the watermarked image. Then, in a step E307, post-processing is applied, such as a 
compression of the digital image. 

[0289] Next, in a step E308, the transformation S is reiterated in order to find the watermarked coefficients, possibly 
made noisy by the post-processing carried out on the image at step E307, and the detectabiiity criterion is calculated 
as previously described. 

[0290] Steps E300 to E309 are reiterated on another acceptable partitioning of the image, here for example by 
dividing the relevant region R into four, and the detectabiiity test is performed at step E309 on the four sub-regions. 
[0291] This partitioning is adopted only if it represents a gain in capacity, that is to say, in this example of quadtree 
division, if the detectabiiity criterion is validated on at least two sub-regions. 

[0292] This variant embodiment makes it possible to directly take into account the post-processings which will be 
applied to the digital image. 

[0293] A method of determining a partition according to a preferred embodiment of the invention will now be de- 
scribed, by reference to Figure 10. 

[0294] The method of inserting a watermark is employed on a digital image I. 

[0295] This insertion method first of all includes a method of determining a partition of a set of coefficients represent- 
ative of the digital image I. 

[0296] This partition is applied as a function of a criterion of detectabiiity of an information bit inserted on each region 
of the partition. 

[0297] It should be noted that, in practice, a partition of a digital image into a maximum number of regions satisfying 
a detectabiiity criterion can be sought so as to ascertain the practical capacity of this image with regard to the insertion 
of a watermark. This maximum capacity of the digital image corresponds to the number of information bits which it is 
possible to insert into the digital image while complying with the detectabiiity criterion. 

[0298] It is also possible, for a given watermark, comprising a given number of information bits of less than the 
maximum capacity of the image I, to seek the best possible partition of the image I which makes it possible to insert 
this watermark of given length, for a given detection probability, and in regions of the largest possible sizes so as best 
to ensure invisibility of the inserted watermark. 

[0299] The method of determining a partition first of all includes a step of spectral decomposition E500 of the wavelet- 
decomposition type described previously (Discrete Wavelet Transform or DWT). 

[0300] Thus a set of spectral coefficients X is obtained, representative of the digital image I. These spectral coeffi- 
cients are distributed into frequency sub-bands as illustrated diagrammatically in Figure 2, on completion of the de- 
composition step E500. The size of the set of coefficients X is equal to N = M X M for a square image. Obviously, the 
same method could be applied to a rectangular image. 

[0301] In this embodiment of the invention, the method includes a step E501 of generating a centred pseudo-random 
sequence of size N equal to the cardinal number N of the set of coefficients X representative of the image I. 
[0302] Let w = {w mn , 0 < m, n <; M} be this pseudo-random sequence. 

[0303] This pseudo-random sequence w is formed from centred pseudo-random sub-sequences w k . 

[0304] The pseudo-random sequence can be seen as the merging of the pseudo-random sub-sequences w k . 

[0305] Here, the centred pseudo-random sequence w is formed from pseudo-random sub-sequences w k of identical 

size which correspond to a block, called unit block, of the set of coefficients X. 

[0306] Obviously, the pseudo-random sub-sequences w k could differ in size from one another. 

[0307] Here, the size of these unit blocks corresponds to the minimum number of coefficients suitable for forming a 

region of a size which is statistically significant for satisfying the detectabiiity criterion over this region. 

[0308] As described later, when a standardised statistical detection test is applied to a set of coefficients, a minimum 
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number of coefficients is necessary for this test to be significant. 

[0309] The pseudo-random sub-sequences w k are of a size A = L x L which is greater than or equal to this minimum 
number cf coefficients. 

[0310] For preference, pseu do- random sub-sequences w*are chosen to be equai in size to this minimum size This 
characteristic is particularly advantageous when partition is achieved by iterative combination by combining blocks 
from an initial segmentation of the image into unit blocks of size L x L. The initial segmentation then corresponds to 
the maximum partition of the image which can be applied, beyond which the blocks of coefficients generated are no 
longer statistically significant for the insertion and detection of an information bit. 
[0311] For example, L = 6 is chosen. 

[0312] The pseudo-random sequence w is then the merging of N/64 pseudo-random sub-sequences w k forminq unit 
blocks. 

[0313] On each unit block, in accordance with the invention, a centred pseudo-random signal of predefined distribu- 
tion (Gaussian, uniform, etc.) has to be generated on the basis of a function depending on a secret key K and on the 
index k of each unit block. This key can be denoted f(K, k). It is then necessary to ensure that, in practice, the average 
of the signal w k is exactly equal to zero over each unit block. 

[0314] A first method for generating such a centred pseudo-random sequence consists in centring the sequence in 
a deterministic way, by symmetrisation. This method is suitable for generating a pseudo-random sub-sequence of even 
size A. 

[0315] In this case, each centred pseudo-random sub-sequence w* is created by generating one half, A/2, of the 
pseudo-random samples of the sub-sequence via a pseudo-random number generator of known centred distribution, 
and by symmetrising the pseudo-random samples thus generated in order to obtain the other half A/2 of the samples 
of the sub-sequence w k . 

[0316] For a single-dimensional sequence, a set of samples of size A/2 {w j( 1 < j < A/2} is therefore generated, and 
this set is made symmetric in order to obtain the other half of the samples {Wj = -w A .j, A/2 < j < A}. 
[0317] In order to reinforce the random character of the sequence thus obtained, permutation can be carried out on 
the set of samples obtained {Wj, 1 < j < A} by using a secret key. 

[031 8] A second method of generating a centred pseudo-random sequence consists in creating each centred pseudo- 
random sub-sequence w k by generating pseudo-random samples via a pseudo-random number generator and by 
redistributing the sum of these pseudo-random samples generated over each sample. 

[0319] Thus, first of all, by using the secret key f(K, k), a sequence of pseudo-random samples of the same size A 

as the unit block, i.e. {w ( -, 1 < i < A}, is generated. 

[0320] Next the exact sum of these samples is calculated: 




[0321] Then this sum S k is redistributed over the samples so as to obtain a centred pseudo-random sequence, that 
is to say one with a sum exactly equal to zero. 
[0322] Thus w*< = (Wj-S^A, 1 < i < A} is obtained. 

[0323] The pseudo-random sequence w is thus composed of the merging of the centred sub-sequences w k in such 

a way that the partitioning of a set of coefficients modulated by the centred pseudo-random sequence w gives rise to 

subsets of coefficients also modulated by centred pseudo-random sub-sequences. 

[0324] Obviously, other techniques could be used to generate centred pseudo-random sequences. 

[0325] As illustrated in Figure 1 0, a step E502 of modulation of the set of coefficients by the centred pseudo-random 

sequence w is then implemented so as to insert an information bit onto the set of coefficients originating from the 

spectral decomposition. 

[0326] In practice, the set of coefficients is copied into a working memory, for example a random-access memory in 
a computer, so as not to mark the image I directly. 

[0327] The same information bit is modulated onto the set of coefficients, for example b = 1 here, according to the 
modulation formula: u 



X'j = Xj + a { Wj, with 0 < i < N. 

[0328] The weighting amplitude Oj is chosen in the usual way so as to guarantee invisibility of the inserted information 
bit. 
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[0329] A constant weighting coefficient can be used for all the coefficients, such that Oj = a v for all i, the value of the 
weighting coefficient a v being equal to a maximum value ensuring imperceptibility of the watermark for modulation of 
the set of coefficients representative of the digital image. 

[0330] It may be worthwhile exploiting the fact that the image signal itself makes it possible to mask the modulation. 
[0331] To this end, it is possible, for each coefficient Xj to be modulated, to use a weighting coefficient according to 
a rule of the form (Xj = kj.cv where kj is a modulation factor depending on the coefficients situated in the vicinity of the 
coefficient in question and cc v is equal to the maximum value of the weighting coefficient ensuring imperceptibility of 
the watermark for the modulation of the set of coefficients. 

[0332] Thus, each coefficient will be modulated as a function of the local content of the image, which makes it possible 
locally to augment the amplitude of the modulation to the benefit of the detectability. 

[0333] For preference, in this embodiment of the invention, it is wished to apply a partitioning to a post-processed 
image signal so as to obtain robustness, a priori of the inserted watermark, against any post-processing undergone 
by the image I. 

[0334] To do that, a distortion is applied, in a distortion step E503, to the set of modulated coefficients. 
[0335] By virtue of the modulation of the coefficients in a single operation, on the basis of a centred pseudo-random 
sequence, the distortion step E503 can be implemented once only upon determination of a partition, in contrast to the 
prior techniques in which each sub-region of coefficients was modulated independently of the other and required dis- 
tortion to be applied to each level of the partition. 

[0336] This is all the more advantageous since, here, the step E503 of applying a distortion includes the following 
sub-steps as illustrated in Figure 10: 

inverse spectral transformation E503a of the set of coefficients X'j modulated in order to obtain a watermarked 
image; 

application proper, E503b, of distortion to the watermarked image; and 

spectral transformation E503c of the watermarked image in order to obtain a set of modulated coefficients after 
distortion. 

[0337] The distortions applied are distortions which alter the value of the pixels without altering the geometry of the 
image, such as the addition of noise, alteration of the contrast, compression, low-pass or high-pass filtering, for exam- 
ple. 

[0338] Compression of the image is considered here, according to the JPEG standard (for Joint Photographic Expert 
Group) with specified quality factor, with q = 75 typically by default for JPEG compression. 

[0339] In practice, an inverse wavelet transformation DWT -1 is applied here, then compression followed by decom- 
pression of the image, and finally forward wavelet transformation again in order to be able to estimate the detectability 
of the information bit inserted in the spectral domain. 

[0340] Next an iterative partitioning of the set of coefficients is applied, so as to search for a partition of the set of 
coefficients in which each inserted information bit can be detected correctly according to a fixed detection probability 
The number of regions thus determined during the partition corresponds to the practical insertion capacity of the image 
I. 

[0341] The means for displaying to the users Q bit determined by calculating the capacity should exist. 
[0342] various techniques for iterative segmentation by division/combination can be used as previously described. 
[0343] A first way of performing the preferred embodiment of the invention will firstly be described, by reference to 
Figure 11, in which the partitioning is achieved by iterative division on the basis of an initial segmentation of the set of 
coefficients into regions of a statistically significant size for validating a detectability criterion over these regions. 
[0344] Following on from the steps of transformation E500, generation E501 , modulation E502 and distortion E503 
as described above, a first frequency sub-band will be considered, in a step E504, as a set of coefficients to be parti- 
tioned. 

[0345] The initial segmentation of this set of coefficients corresponds to a spatial segmentation of the image I into 
blocks of unit size as defined above. 

[0346] By its very construction, the pseudo-random modulation sequence is centred on each of these unit blocks. 
[0347] In a quadtree-mode combination, it is always four blocks of equal size which are grouped together, called 
"child blocks' in order to form a "parent block", and this combination is recommenced iteratively giving the notion of 
recursiveness. 

[0348] Obviously, partition takes place frequency sub-band by frequency sub-band of the signal arising from the 
spectral decomposition. 

[0349] Each unit block is considered here, in a step E505, in an order of scanning of the sub-band, for example from 
left to right and from bottom to top. 

[0350] A validity criterion is applied over a first unit block which has not yet been processed. 
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[0351] In a calculation step E506, a statistical detection test is applied to this unit block of modulated coefficients XV 
[0352] An example of a standardised statistical test is given below. 

[0353] Let X\, 1 < i < P be the set of coefficients to be tested after insertion of the pseudo-random sequence w and 
application of pre -distortion. 

[0354] The test consists in calculating, for this set of size P: 

T-™JP 

in which M represents an average value of correlation between the values of each coefficient X'j and of modulation 
w h and V the variance of these values. 
[0355] Hence in practice: 

M = i f w. (x f . —M x ) , With M x = UH± 

P i~\ * p 

and 



[0356] The calculation of this test value T requires knowledge of the pseudo-random modulation sequence w which 
is generated once and for all at the above-described generation step E501 . 

[0357] The detectability is given by the comparison of this test value T with a certain threshold P which corresponds 
theoretically to a probability of correct detection of p% if the test T follows a Gaussian law with variance equal to 1 for 
the two hypotheses tested (presence or absence of watermark). 

[0358] A comparison step E507 makes it possible to compare this test value T with the threshold value P or correct 
detection rate. p 

[0359] On completion of this comparison, if the calculated value T is higher than the correct detection rate, the block 
tested is held in memory, in a memory-storage step E509, as being a region available for inserting an information bit. 
[0360] In practice, a capacity-calculating step E510 can be implemented in such a way that, for every block of coef- 
ficients tested which is suitable for having an information bit inserted, the value of the capacity is increased by one unit. 
[0361] In contrast, if the calculated value T is less than the correct detection rate P p , a block of higher level is selected, 
in a combining step E50B, that is to say, here, a parent block corresponding to the combination of four unit blocks, and 
the steps E506 of calculating the test value T and of comparison E507 with the correct detection rate P D are reiterated 
over this new block. 

[0362] During this combining of child blocks in order to form a parent block, those child blocks are preferably chosen 
for which the detectability criterion is not satisfied. 

[0363] For application of this detectability criterion, any standardised statistical test can be used for the detection, 
such as that described in the article entitled "A method for signature casting on digital images" by I. PITAS, in Proc.' 
ICIR pages 215-218, September 1 996, which makes it possible to characterise the detection in terms of probability 
[0364] A threshold value P p can be chosen, corresponding to a fixed level of detection probability, 99.95% for ex- 
ample. 

[0365] Hence, for a given region defining a set of coefficients X of size P and with a fixed weighting coefficient a, it 
is possible to calculate a detectability amplitude and to compare the latter with a threshold value P so as to validate 
or not validate the detectability criterion over this region. P 
[0366] In a general way, the detection probability decreases with the size of the regions tested, such that, by com- 
bining regions in order to increase the size of the block of coefficients used for inserting a watermark, the value of the 
test T is increased over this new block. 

[0367] Hence, on completion of these steps of adaptive partitioning E504 to E508, an image insertion capacity can 
be determined, equal to Q bits, as can an appropriate partition of the image corresponding to this capacity Q. 
[0368] A watermark can then be formed, with length equal to Q bits, and can then be inserted, bit by bit, into the 



20 



EP 1 043 687 A1 

various regions of the partition. 

[0369] Obviously, the step E520 of insertion proper of the watermark can be coupled with the determination of the 
partition, in such a way that each region available on conclusion of the detectability test E506 E507 can actually be 
watermarked by the insertion of an information bit by modulation of the coefficients of this region 
[0370] The msertion step ES20 makes it possible actually to insert an information bit of the watermark by choosina 
a coefficient value b = 1 or = -1 in the above^iescribed modulation formula, the alternating of the bits conveying 
the content of the message inserted. a 
[0371] A centred pseudo-random sequence of known distribution can be used, different for each region of the oar- 
tition, and initialised by a secret key K', as illustrated in Figure 1 . 

[0372] It should be noted that this secret key K' may, if appropriate, be identical to the secret key K used during the 
generation of the pseudo-random sequence w for determining the partition. 

[0373] Moreover, a centred pseudo-random sequence w identical to that generated at the generation step E501 
could be generated. K 

L 0 ^ 7 ?] «n= rtili0 u 0 L a ? im39e ' iS Sh ° Wn dia 9rammatically by way of example in Figure 10. on conclusion of steps 
f™° ~ " example ' 25 se Parate regions are obtained making it possible to insert a 25-bit watermark 

10375] Obviously, any other partitioning technique could be used. 

[0376] A second way of performing the preferred embodiment of the invention will now be described by reference 
to Figure 1 2, in which the partitioning is achieved by iterative division of the set of coefficients 

[0377] On conclusion of the steps of transformation E500, generation E501. modulation E502 and distortion E503 
fn,^ " * „° V6, 3 ,requenc y sub " ban d will be considered, in a step E51 1 . as set of coefficients to be partitioned 
[0378] A quadtree-mode division will be considered here, in which one block, called 'parent block" is always divided 
into four blocks of equal size, called "child blocks". 

[0379] Obviously, partition takes place frequency sub-band by frequency sub-band of signal originating from the 
spectral decomposition. 

[0380] Hence, in a step E51 2. an initial block corresponding to a first frequency sub-band is considered 
[0381] This initial block is stored in memory in a list L of blocks to be processed 

conSered 6 *' * & bl °° k B ^ ^ prOCeSSSd ' COmin9 from ,he list °' blocks to be processed L. is 

[0383] A detectability criterion is applied on this block B in a checking step E514 

[0384] In practice, as previously described, the value T of a statistical detection test on this block of modulated 
coefficients is calculated in a calculation step E514a. This step is analogous to the calculation step E506 described by 
reference to Figure 11. ' 

[0385] The detectability is then given by the comparison of this test value T with a certain threshold P B . which theo- 
retically corresponds to a correct detection probability of p% if the test T follows a Gaussian law with variance equal 
to 1 for the two hypotheses tested (presence or absence of watermark). 

[0386] A comparison step E514b makes it possible to compare this test value T with the threshold value P or correct 
detection rate. This comparison step E51 4b is analogous to the comparison step E507 described by reference to Figure 

[0387] On completion of this comparison, if the calculated value T is less than the correct detection rate P a check 
is made, in a test step E515, as to whether there are blocks remaining to be processed in the list L 
[0388] If not, the recursive partitioning algorithm is terminated. 

[0389] Otherwise the following block B of the list L is considered at step E513, and the checking step E51 4 is reit- 

013100. 

[0390] If, on completion of the comparison step E514b, the calculated value of the test T is higher than the correct 
detection rate P p . detectability is tested over each child block of this parent block B in a checking step E51 6 
[0391] In practice, the parent block B is divided into four child blocks and the detectability criterion is checked on 
each of these blocks according to a step analogous to the previously described checking step E514 

A ' eSt E51 7 mak6S ' l possib,e to verif y wh ether the detectability criterion is valid over at least one child block 
[0393] If yes, a division step E518 is implemented in order to actually divide the block B into four child blocks and 
an addition step E51 9 allows these child blocks to be added to the list of blocks to be processed L 
[0394] The test and checking steps E51 3 to E51 9 are then reiterated over each of the blocks of the list L until the 
latter is empty. 

[0395] In contrast, if. on conclusion of test E51 7. the detectability criterion is not satisfied on any of the child blocks 
of the parent block B. the tested block B is held in memoiy in a memory-storage step E509 as a region available for 
inserting an information bit. " 

[0396] In practice, a step of updating the capacity E510 can be implemented as described previously, in such a way 
that, for each block of coefficients tested which is suitable for insertion of an information bit. the value of the capacity 
is increased by one unit. 



21 



EP 1 043 687 A1 



[0397] Hence, on conclusion of these adaptive partitioning steps E511 to E519, an image insertion capacity can be 
determined, equal to Q bits, as can an appropriate partition of the image corresponding to this capacity Q. 
[0398] A watermark can then be formed, with length equal to Q bits, and can then be inserted, bit by bit, into the 
various regions of the partition. 

[0399] Obviously, as before, the step E520 of insertion proper of the watermark can be coupled with the determination 
of the partition, in such a way that each block B available on conclusion of the detectability test E507, can actually be 
watermarked by the insertion of an information bit by modulation of the coefficients cf this block B. 
[0400] The insertion step E520 is carried out in the same way as when partition is achieved by iterative combination. 
[0401] Moreover, in one enhanced embodiment of this technique of partition by iterative division, a division of a parent 
block B into child blocks can be adopted on completion of the detectability test E507 only if it represents a gain in 
capacity, that is to say into a quadtree-mode division only if the detectability criterion is satisfied over at least two child 
blocks. Otherwise, it may be preferable to keep the parent block of larger size so as to guarantee better invisibility of 
the inserted watermark. 

[0402] Obviously, the partitioning techniques described previously are not in any way limiting. 
[0403] Notably, it is possible to use various iterative segmentation techniques by division/merging, to which there 
belong the quadtree block partitioning described previously, but also graph partitioning. In general terms, there are the 
so-called bottom-up methods and the top-down methods depending on whether the iterative partitioning means is 
merging or division. 

[0404] Hence it is possible, if appropriate, to combine two techniques of partitioning by division and combination. 
This is the case when an initial partition is chosen, for example, intermediate between the two extremes which are, on 
the one hand, a single region and, on the other hand, a partition into regions of the smallest admissible size, allowing 
both the operations of division and of combination. 

[0405] It is also possible to implement the method of determining a partition in order to maximise the capacity of the 
digital image by using several different segmentation techniques and then choosing the partition which confers the 
greatest capacity Q for the digital image. 

[0406] It is also possible to determine the best partition possible, that is to say the one yielding a set of blocks of 
coefficients of largest possible size, for inserting a fixed-length message, less than the maximum capacity of the digital 
image to be watermarked. 

[0407] In all cases, the present invention makes it possible to generate and to insert a pseudo-random sequence 
once only over a set of coefficients which are available for carrying out the insertion, and, as a corollary, to apply only 
one single pre -distort ion to the set of coefficients. 

[0408] It is thus possible rapidly to obtain a partition of the set of coefficients which makes it possible to insert a 
watermark with a fixed correct detection probability for an entire range of distortion lying between "no distortion" and 
•predetermined distortion" applied during the determination of the partition. Moreover, the simulation of maximum dis- 
tortion before insertion is a quasi-deterministic operation which confers an advantage with respect to a theoretical 
statistical calculation. 

[0409] A description will be given now of the insertion method according to the invention in a preferred embodiment. 
[0410] The watermarking signal is inserted after a transformation step E600, in a set of coefficients of the domain 
transformed into sub-bands. 

[0411] As illustrated in Figure 13, in a choosing step E601, there is for example chosen the high-frequency sub-band 
HH V with the first resolution level, corresponding to a high-pass filtering in the horizontal and vertical directions. 
[0412] There is thus a set of spatio-frequency coefficients, with a cardinal number equal to N, denoted for example 
X =< Xj, 1 < i < Nr. 

[0413] For an image I of size 512X512, the sub-band HH 1 is of size N = 256X256. 

[0414] Since the watermarking has to be imperceptible and indelible, and therefore difficult to locate and remove by 
pirating, it is chosen to insert a pseudo-random watermarking signal, spreading its spectrum in order to make this signal 
undetectable by spectral or statistical analysis. 

[0415] For example, a pseudo-random sequence w is considered, which follows a uniform law on the interval [-1,1] 
with w -\ w jt 1 < i < P>, the length P of the sequence being less than or equal to N. 

[0416] Naturally, any pseudo-random watermarking signal, with a known distribution and a nil mean, may be suitable. 
The most usual distributions for the watermarking signal w are, apart from the uniform distribution on [-1 , 1] mentioned 
above, the binary distribution <-1 ,1r and the centred standardised Gaussian distribution N(0, 1). 
[0417] The modulation which it is wished to apply to the coefficients in order to insert the watermarking signal uses 
a linear model of the type: 

X i = X j(«) + aw i' with 1 - j - p where 
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[0418] X j(j) is a subset of spectral coefficients chosen from the set of coefficients X, and a is a weighting coefficient. 
[0419] Here, the weighting coefficient a is constant for any j. 

[0420] In this example embodiment, it is sought to determine whether the insertion of the watermarking signal w is 
actually possible in the set of coefficients X whilst meeting predetermined criteria of imperceptibility and detectability. 
5 [0421] According to the invention, the insertion method first of all includes a step E602 of calculating a maximum 
value oty of the weighting coefficient ocas a function of the length P of the watermarking signal w ensuring imperceptibility 
of the watermarking signal w. 

[0422] This maximum value or visual amplitude corresponds to the maximum modulation amplitude which can 
be used in the aforementioned linear insertion model beyond which an observer is capable of visually detecting a 

io change at the level of the reconstructed watermarked image. 

[0423] In this embodiment, a visibility model is used which makes it possible to predict the visibility of a watermarking 
operation as a function of the different parameters, which are the representation of the signal through the spatio- 
frequency transformation T used, the sub-band HH-, considered for the insertion, the type of distribution of the sequence 
of the watermarking signal w and the length P of the sequence w 

15 [0424] A simple model, developed by WATSON and described in the article entitled "Visibility of wavelet quantization 
noise", A.B. WATSON et al, IEEE Trans, on Image Process, 6(8), 1164-1175, 1997, makes it possible to predict the 
visibility of a set of coefficients modulated from the measurement of visibility of a single modulated coefficient. Reference 
can advantageously be made to this document for a detailed description of this model. 

[0425] During the step E602 of calculating the maximum value oty, use is thus made of a function dependent on the 
20 length P of the watermarking signal w and on the type of spectral transformation T used, but independent of the set of 
coefficients to be watermarked X. 
[0426] This calculation function used can be written: 



25 a v (P,S,w)= ^ ' 

p *(E[\wr]y p 

where 

30 ot tase( s ) is a Dase value, dependent on the transformation used S and the (base) sub-band considered for the 

insertion, of the maximum weighting coefficient ensuring imperceptibility when a single coefficient of this sub-band 
is modulated, 

P is strictly greater than 2, and 

£[lwfp is the mathematical expectation of the function IwlP. 
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[0427] The base values 04,^(5) can be measured once and for all, for each sub-band of coefficients in the wavelet 

decomposition, from a single psychovisual measurement, and be stored in a table of visual amplitudes. 

[0428] p is the exponent of a Minkowsky sum and can be chosen equal to 5, for example. 

[0429] The mathematical expectation EflwP] corresponds to an estimation of the mean of the function IwlP. 

[0430] This visibility model does not take into account the image I to be watermarked itself and is independent thereof. 

This is equivalent to considering that the image I is uniform. It is a question of a "worst case" model since the presence 

of the image signal makes it possible to visually mask the modulation itself. 

[0431] In this example, this calculation step E602 is implemented in order to calculate the maximum value of the 
weighting coefficient a for a length P of the watermarking signal w equal to the cardinal number N of the set of modu- 
latable coefficients X. Thus a maximum value denoted a^N) is calculated, which ensures imperceptibility of the wa- 
termarking signal w at the limit of perceptibility. 

[0432] In the example embodiment described here, for a value of p = 5, the visual amplitude for the set of coefficients 
of the sub-band HH, is given by: 



50 



a v (N) = 



6 a 



HH1 



(7) 



[0433] Next, in a calculation step E603, the length P of the watermarking signal w is calculated for a minimum value 
Oq of the weighting coefficient a equal to said maximum value a^N) calculated and for a predetermined detection 
probability level. 

[0434] This is because, in order to meet a predetermined detectability criterion, the weighting coefficient a must be 
greater than a minimum value a D referred to as the detection amplitude. 
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[0435] At the decoder, the received image r, which corresponds to the watermarked initial image, possibly made 
noisy because of certain distortions applied to the image, is first of ali transformed by a transformation identical to that 
used during the insertion of the watermarking signal in order to find once again the set of coefficients which are mod- 
ulated. 

[0436] Thus a set X* =< X^, 1 < i < N) is obtained, which represents the set of coefficients liable to have been modulated 
By hypothesis, it is assumed that any distortion suffered by these coefficients can be modelled by an additive noise 
decorrelated from the image signal itself n ={n r 1 < i < Nr.. 
[0437] Each coefficient can thus be written: 

Xj* = Xj + awj + n jt with 1 < i < N. 

[0438] The detection of the inserted watermarking signal w consists of constructing a hypothesis test for replying to 
the question "has the signal w been inserted in X*?". 

[0439] The advantage of constructing a hypothesis test is that its result makes it possible to give a theoretical con- 
fidence level to the decision since a hypothesis test is standardised with respect to known distributions. 
[0440] For a class of statistical tests based on correlation, it can be shown that this detection amplitude can be 
calculated from a calculation function of the type: 



I 2 2 

a(Jo y + oj 

where a* is equal to the variance of the coefficients to be modulated, o* is equal to the variance of the additive 
noise and a, b and c are constants which depend on the distribution of the watermarking signal w and on the required 
detection probability level p. 

[0441] Thus this minimum value Oq is determined according notably to the length P of the watermarking signal w. a 
detection probability level p, the variance of the coefficients to be modulated, the distribution of the watermarking signal 
and possibly the variance of an additive noise decorrelated from the coefficients modelling any distortions applied to 
the coefficients. 

[0442] Conversely, it is also possible to express the minimum detection length, denoted P°, as a function of the 
weighting coefficient a, for a fixed detection probability level: 

2 2 2 2 

da 

[0443] It would be possible to use this result to check the parameters which make it possible to insert a watermarking 
signal w which will be robust to any additive noise with an energy less than or equal to a n 2 . 

[0444] As illustrated on the curves in Figure 14, from the visual amplitude a v calculated as a function of the length 
P of the watermarking signal to be inserted, the insertion method includes, in general terms, a step of determining 
compatible values of the length P of the watermarking signal w, of the weighting coefficient a and of the probability 
level p so that the value of the weighting coefficient determined a is less than or equal to the maximum value or visual 
amplitude cc v and greater than or equal to the minimum value or detection amplitude for these determined values 
of the length P of the watermarking signal w and the detection probability level p. 

[0445] The weighting coefficient a, the length P of the watermarking signal w and the detection probability level p 
must thus be determined so that the following inequality is satisfied: 



[0446] As illustrated in Figure 14, in order to calculate the minimum value P° of the watermarking signal w for which 
the aforementioned inequality will be satisfied, the formula linking the length P of the watermarking signal w to the 
amplitude of the weighting coefficient a is used, with a = a,/N) and with a detection probability level of p% 
[0447] In this example a standardised statistical test is taken, based on the correlation between the watermarked 
signal X* and the inserted watermarking signal w, as proposed for example in the document entitled "On Resolving 
Rightful Ownerships of Digital Images by Invisible Watermarks" by W. ZENG and B LIU in Proc ICIP 97 D aaes 
552-555, October 1997. 
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[0448] It can be shown that, for this test: 



[0449] A calculation step E604 thus makes it possible to calculate the three aforementioned constants a b and c so 
that 



a=t 0 JE( W 2 ), 



and 



c= $o 2 2 . 

w 

[0450] The values E(vfi)ar\6o 2 2 can be calculated theoretically as a function of the distribution of w. For the uniform 
distribution considered here, the Theoretical values are E(w*)= 1/3 and o 2 2 = 4/45. 

[0451] The value ^ is a threshold value, dependent on the probability oT correct detection p : t 0 = 2 x t . The values 
of t p are percentile of a centred Gaussian law and have a variance equal to 1 . These values are tabulated in the state 
of the art, for example in the book entitled "Probability and Statistics" by A. PAPOULIS, Prentice-Hall, 1990. 
[0452] For example, t p = 3,291 corresponds to p = 0,9995, that is to say to a detection probability level equal to 
99,95%. 

[0453] In this numerical example it is considered that any additive noise is not taken into account so that o n = 0. 
[0454] In an estimation step E605, the variance of the set of coefficients X is also estimated. The variance is estimated 
according to the following formula: 



with 



1 N 

[0455] Then, for a detection probability fixed at 99.95%, the function linking the length P of the watermarking signal 
w to be inserted and the weighting coefficient a is then derived therefrom: 

14,44 X ol+ 3.84a 2 



[0456] At the calculation step E603, the minimum length P° is determined from the above formula for a weighting 
coefficient equal to the maximum value <x v (N) calculated at the calculation step E602. 

[0457] The insertion method includes, in this example embodiment, a step E606 of comparing this calculated length 
P° with the cardinal number N of the set of modulatable coefficients X. 

[0458] If, at this comparison step E606, the calculated length P° is less than or equal to the cardinal number N of 
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the set of modulatable coefficients X, it is deduced therefrom that the insertion of a watermarking signal w of length P> 
is possible and meets the required imperceptibility and detection criteria. 

[0459] In reality, it is possible to effect an insertion of a watermarking signal w which meets these criteria since the 
length P of this watermarking signal lies between the calculated value P° and the cardinal number N of the set of 
modulatable coefficients X. 

[0460] Then the insertion proper of a watermarking signal of length P is carried out by modulating a subset of coef- 
ficients with a cardinal number P according to the linear modulation level, with a weighting coefficient a equal to the 
calculated maximum value a v (N). 

[0461] The insertion method includes, in the usual fashion, a reverse spatio-frequency transformation step which 
makes it possible to obtain, at the output of the coder, the watermarked image I*. 

[0462] Where insertion is possible, it is also possible to calculate the number of theoretical bits which can be inserted 
This is because, at the decoder, the decoding operation consists of detecting whether a pseudo-random sequence w 
has been inserted. The response being binary (yes/no), it can be considered that a pseudo-random sequence of length 
P makes it possible to insert and extract an information bit in the image. By repeating the insertion operation on several 
subsets of spectral coefficients, it is thus possible to insert and extract several information bits. 
[0463] Consequently, since in this example there are a set of coefficients of cardinal number N and since only a 
number P° of coefficients are necessary to insert and detect an information bit, it is possible to derive therefrom the 
total number Q of information bits which can be inserted, referred to as the capacity of the image* Q = N/P° 
[0464] Moreover, if at the comparison step E606 the calculated length P° is greater than the cardinal number N of 
the set of modulatable coefficients X, it is a priori deduced therefrom that the insertion of a watermarking signal meeting 
the criteria of imperceptibility and detectability is not possible. 

[0465] In this case, the insertion method can if necessary include an additional step of calculating the detection 
probability level p' obtained for a minimum value a D of the weighting coefficient equal to the maximum value a v (N) 
calculated for a length P of the watermarking signal w equal to the cardinal number N of the set of modulatable coef- 
ficients X. 

[0466] If this detection probability level, which will be less than that fixed initially at 99.95%, is judged satisfactory, 
the insertion of a watermarking signal of length equal to N can be effected, that is to say all the coefficients of the sub- 
band of frequency HH, are modulated. 

[0467] In a second embodiment of the invention, the insertion method can include a step of calculating a threshold 
value P* of the length of the watermarking signal w, as illustrated in Figure 14. 

[0468] This threshold value P* is determined so that the minimum value a 0 of the weighting coefficient is equal to 
the max.mum value a v of the weighting coefficient for a length P of the watermarking signal w equal to this threshold 
value P* and a predetermined detection probability level p. 

[0469] This calculation step makes it possible to determine, from formulae expressing the detection amplitude a D 
and the visual amplitude cc v notably as a function of the length P of the watermarking signal w, the minimum value P* 
of the insertion signal from which the weighting coefficient a can satisfy the inequality ct D < a <a v . 
[0470] In Figure 14, this minimum value P* is given by the intersection of the two curves representing the variation 
in the visual amplitude a v and the variation in the detection amplitude as a function of the length P of the inserted 
watermarking signal w, for a predetermined detection probability level. 

[0471] Thus, in general terms, the insertion method according to the invention has the advantage of being able to 
control a priori and theoretically the detection probability and imperceptibility of the watermarking signal which it is 
sought to insert in an image. 

[0472] Naturally, many modifications could be made to the example embodiments described above without departinq 
from the scope of the invention. 

[0473] Thus the insertion technique used could be applied to the raw digital image, without undergoing any spatio- 
frequency transformation prior to the modulation of the coefficients. 

[0474] In this case, the modulated coefficients are coefficients representing the digital image solely in the spatial 
domain. 

[0475] Moreover, the spatio-frequency transformation applied to the image can use analysis and synthesis filters 
other than those described above, or even be a transformation other than discrete Fourier transformation discrete 
cosine transformation or Fourier-Mellin transformation. These transformations are used currently in conventional 
processings of digital images. 

[0476] In addition, the invention is not limited to the example embodiments described but also concerns any embod- 
iment where the values of the weighting coefficient a, of the length P of the watermarking signal and of the detection 
probability level p are determined in a correlated fashion so as to satisfy the inequality oq £ a <cv 
[0477] The present invention also concerns a device for determining a partitioning and a decoding device which are 
adapted to implement the partitioning and decoding methods described previously. 

[0478] These devices for determining a partitioning and decoding can be used in a computer 50 as illustrated in 
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Figure 15, independently or in the same computer 50. 

[0479] The partitioning determination device comprises partitioning means 500, 502, 503 adapted to perform an 
adaptive partitioning according to a criterion ot detectability of an information bit inserted on each region. 
[0480] These partitioning means use on coding a maximum value <x v of the weighting coefficient ensuring impercep- 
tibility of the watermarking signal for modulating all the coefficients representing said digital image, possibly with a 
modulation factor kj dependent on each coefficient, as described previously. 

[0481] This device for determining a partitioning also has means 500, 502, 503 of comparing the size P of each 
region of the partitioning with a minimum size P stat corresponding to minimum size of a statistically significant sample 
for the detection of a watermarking signal inserted in said region. 

[0482] The device also has means 500, 502, 503 for applying a pre-distortion before calculating a detectability cri- 
terion. 

[0483] The partitioning means comprise means 500, 502, 503 of modulating all the coefficients of each region by 
inserting the watermarking signal, means 500, 502, 503 of calculating a detectability amplitude from the detectability 
criterion and means 500, 502, 503 of validating the detectability criterion, adapted to compare the detectability amplitude 
with a predetermined threshold value Tc. 

[0484] The partitioning means 500, 502, 503 are adapted to effect a partitioning by iterative division of the numerical 
image and, for each region at a rank in the said partitioning, the validation means 500, 502, 503 are adapted to validate 
a partitioning of a rank higher in the region if and only if there exist at least two subregions of the region for which the 
detectability criterion is validated. 

[0485] This device also comprises means 500, 502, 503 of incrementing the value of a capacity Q of the digital image, 
adapted, for each region at a rank in the partitioning, to increment the capacity value Q when the partitioning of rank 
directly higher is not validated, the capacity Q being equal to the cardinal number of the distinct regions obtained by 
the partitioning on which the detectability criterion is validated. 

[0486] In another embodiment, this device for determining a partitioning comprises incrementation means 500, 502, 
503 adapted, for each region R at a rank in the partitioning, to increment the value of a capacity Q when the partitioning 
of rank directly higher is not validated, comparison means 500, 502, 503 adapted to compare the value of the capacity 
Q with a fixed capacity value Qs and reiteration means 500, 502, 503 adapted to reiterate the partitioning for a region 
having the highest detectability amongst all the other regions to be processed when the capacity value Q is less than 
the fixed capacity value Qs. 

[0487] In a third embodiment of the invention, the partitioning means 500, 502, 503 are adapted to effect a partitioning 
by iterative merging of the digital image and there are adapted, for each region of the partitioning for which the de- 
tectability criterion is not validated, to merge the region with another region of the partitioning, and preferably with 
another region of the partitioning for which the detectability criterion is not validated. 

[0488] In an entirely similar fashion, the device for decoding a watermarking signal in a digital image comprises 
means 500, 502, 503 of partitioning the digital image to be decoded into distinct regions, adapted to perform an adaptive 
partitioning as a function of a criterion of detectability of an information bit obtained by demodulating the coefficients 
on each region. 

[0489] it also has means 500, 502, 503 of comparing the size P of each region R of the partitioning with a minimum 
si2e p stat corresponding to the minimum size of a statistically significant sample for the detection of a watermarking 
signal inserted in the region. 

[0490] It comprises means 500, 502, 503 of calculating a detection amplitude from the detectability criterion and 
means 500, 502, 503 of validating the detectability criterion, adapted to compare the detection amplitude with a pre- 
determined threshold value Td for decoding. 

[0491] This predetermined threshold value Td for decoding is less than the predetermined threshold value Tc used 
during the partitioning determination method described previously. 

[0492] In a preferred embodiment, the device for determining a partition includes means 500, 502, 503 for checking 
a criterion of detectability, over each region, of an information bit inserted by means of a pseudo-random sequence 
with average of zero by modulation of the coefficients of the region. 

[0493] It also includes means 500, 502, 503 for generating a centred pseudo-random sequence equal in size to the 
cardinal number of the set of coefficients representative of an image, here equal to N. 

[0494] Means 500, 502, 503 of modulating the set of coefficients X by the centred pseudo-random sequence in order 
to insert the same information bit onto this set of coefficients are also incorporated into the device for determining a 
partition 11 as shown on figure 1. 

[0495] The latter may also includes means 500, 502, 503 for applying a distortion D which makes it possible to 
provide for and guarantee the robustness of the inserted watermark in the case of certain types of distortion, typically 
a compression of the digital image in order to store it. 

[0496] These means of applying a distortion make it possible actually to simulate a predetermined distortion over 
the set of modulated coefficients. Here they comprise, whenever the coefficients of the set X are spectral coefficients, 
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means 500, 502, 503 for inverse spectral transformation of the set of coefficients modulated in order to obtain a wa- 
termarked image, means 500, 502, 503 for applying a distortion to this watermarked image, and means 500, 502 503 
for spectral transformation of the watermarked image in order to recover the set of modulated coefficients after distor- 
tion. 

[0497] In this example, the inverse -spectral-transformation means are suitable for implementing multi-resolution 
spectral wavelet recomposition and the forwara spectral transformation means are suitable for again reiterating the 
breakdown of the digital image I into wavelets. 

[0498] The device for determining a partition 11 , in a first way, also includes means for partitioning by iterative com- 
bination, on the basis of an initial segmentation of the set of coefficients X. 

[0499] For preference, this initial segmentation makes it possible to segment the set of coefficients X into a set of 

regions of a size which is statistically significant for validating the delegability criterion over these regions 

[0500] The means of partitioning 500, 502, 503 by iterative combination are suitable, for each region of the partition 

for which the detectability criterion is not satisfied, for combining the region with another region of the partition 

[0501] Obviously, the partitioning means 500, 502, 503 could also be suitable for forming this partition of the image 

by iterative division of the digital image according to a second way of performing the preferred embodiment of the 

invention. 

[0502] This partition-determining device 1 1 makes it possible, as output, to obtain a partition into regions of the set 
of coefficients X, with it being possible to insert an information bit onto each of these regions with a reliable detection 
probability and robustness toa given distortion. The maximum number of regions which can be obtained by this adaptive 
partitioning of the set X corresponds to the practical insertion capacity of this set X. 

[0503] As illustrated in Figure 1 , an insertion device 10, incorporating the device for determining a partition 11 de- 
scribed above, moreover comprises means 1 3 for reading a watermark to be inserted, containing several information 
bits. By way of an example which is not at all limiting, a binary watermark may be 1000111 . representing, for example, 
the coding of the name of the author of the image to be watermarked. 

[0504] The insertion means 1 2 proper may be conventional and make it possible, by modulation, to insert the various 
information bits of the watermark onto the various insertion carriers determined by the device for determining a partition 

[0505] According to the invention, the device for inserting a watermarking signal w comprises means 500, 502, 503 
of spatio-frequency transformation of an image I, and for example analysis filters associated with decimators by two 
adapted to effect a wavelet decomposition of an image I. It also has reverse spatio-frequency recomposition means 
500, 502, 503 for recomposing the image I after the insertion of the watermarking signal in the domain transformed 
into sub-bands. 

[0506] It also comprises means 500, 502, 503 of inserting a watermarking signal w of length P adapted to modulate 
a subset of coefficients of cardinal number P according to the linear model described previously using a weiqhtinq 
coefficient a. a a 

[0507] This insertion device also comprises means 500, 502, 503 of calculating the maximum value a v of the weight- 
ing coefficient as a function of the length P of the watermarking signal w ensuring imperceptibility of the watermarking 
signal. These calculation means use the function a v (P, T) described above and dependent on the length P of the 
watermarking signal w and the spectral transformation T used for decomposing the image signal I. 
[0508] In general terms, this insertion device also has means 500, 502, 503 of determining compatible values of the 
length P of the watermarking signal w, the weighting coefficient a and the probability level p so that the value of the 
weighting coefficient determined is less than or equal to the maximum value <x v and greater than or equal to the minimum 
value ceo for the determined values of the length P of the watermarking signal w and the detection probability level p. 
[0509] The calculation means 500, 502, 503 are also adapted to determine a minimum value or detection amplitude 
otp using the calculation functioning described above of the type: 



2 2 . 
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[0510] In one embodiment of the invention, the insertion device also has means of calculating 500, 502, 503 the 
length P° of the watermarking signal w for a minimum value cc D of the weighting coefficient equal to a maximum value 
cc v (N) calculated for a length P of the watermarking signals w equal to the cardinal number N of the set of modulatable 
coefficients ensuring imperceptibility of the watermarking signal w at the limit of perceptibility and for a predetermined 
detection probability level p; and means 500, 502, 503 of comparing this length P° calculated with the cardinal number 
N of the set of modulatable coefficients X. 

[0511] It also comprises means 500, 502, 503 of calculating a detection probability level p' obtained for a minimum 
value Oq of the weighting coefficient equal to the maximum value a v calculated for a length P of the watermarking 
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signal equal to the cardinal number N of the set of modulatable coefficients X. 

[0512] In another embodiment, the insertion device comprises means 500, 502, 503 of calculating a threshold value 
P* of the length of the watermarking signal w determined so that the minimum value cc D of the weighting coefficient is 
equal to the maximum value a v of the weighting coefficient of the length P of the watermarking signal w equal to this 
threshold value P* and a predetermined detection probability level p. 

[0513] All the previously described devices are incorporated in a microprocessor 500 of the computer 50, a read- 
only memory 502 being adapted to store a program for determining a partitioning and/or inserting a watermarking 
signal and/or decoding the watermarking signal and a random access memory 503 comprising registers adapted to 
store variables modified during the running of the programs. 

[0514] The microprocessor 500 is integrated into the computer 50, which can be connected to different peripherals 
for example a digital camera 507 or a microphone 511 , by means of an input/output card 506 in order to receive and 
store documents. 

[0515] The digital camera 507 makes it possible notably to supply images to be authenticated by insertinq a water- 
marking signal. 

[0516] This computer 50 has a communication interface 512 connected to a communication network 51 3 in order to 
receive, if necessary, images to be watermarked. 

[0517] The computer 50 also has document storage means, such as a hard disk 508, or is adapted to cooperate by 
means of a disk drive 509, with removable document storage means such as diskettes 510. 

[0518] These fixed or removable storage means can also contain the code of the different methods according to the 
invention which, once read by the microprocessor 500, will be stored on the hard disk 508. 

[0519] By way of variant, the program enabling the devices to implement the invention can be stored in the read onlv 
memory 502. 

[0520] As a second variant, the program can be received in order to be stored as described above by means of the 
communication network 513. 

[0521] The computer 50 also has a screen 504 for serving, for example, as an interface with an operator by means 
of the keyboard 514 or any other means. 

[0522] The central unit 500 (CPU) will execute the instructions relating to the implementation of the invention. On 
powering up, the programs and methods relating to the invention stored with a non-volatile memory, for example the 
read only memory 502, are transferred into the random access memory 503 (RAM) which will then contain the exe- 
cutable code of the invention and the variables necessary for implementing the invention. 

[0523] This random access memory 503 contains different registers for storing the variables necessary for the running 
of the program, and notably a register for storing the coefficients of the regions R which are temporarily modulated or 
demodulated in order to verify the detectability criterion, the minimum value P stat , and the threshold value Tc and Td 
respectively on coding and decoding. 

[0524] This random-access memory 503 also includes an X i( Xj register for storing the coefficients of the regions 
which are modulated or demodulated temporarily so as to check the detectability criterion and a register w for storing 
the centred pseudo-random sequence generated in order to determine the partition of the set X. 
[0525] Finally, the random access memory 503 contains a register for storing the table of visual amplitudes of 
(T), a register for storing the spectral coefficients X, a register for calculating the visual amplitude a v (N) at the limitof 
imperceptibility, a register for storing the values of the constants a, b, c and of the variance cy* necessary for determining 
the detection amplitude oc D as a function of the length P of the watermarking signal w and a register for storinq the 
minimum value P° calculated. 

[0526] A communication bus 501 affords communication between the different sub-elements of the computer 50 or 
linked to it. The representation of the bus 501 is not limitative and notably the microprocessor 500 is able to indicate 
instructions to any sub-element directly or by means of another sub-element. 



Claims 



2. 



Method of determining a partitioning into distinct regions of a digital image (I) in order to insert a watermarking 
signal (W) by modulating coefficients (X) of said regions (R) representing a digital image (I), characterised in that 
the partitioning into regions is effected by an adaptive partitioning as a function of a criterion of detectability of an 
information bit inserted on each region (R). 

Method of determining a partitioning according to Claim 1 , characterised in that the detectability criterion is calcu- 
lated by using a maximum value (a v ) of the weighting coefficient ensuring imperceptibility of the watermarking 
signal for modulating all the coefficients representing said digital image. 
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3. Method of determining a partitioning according to claim 2, characterised in that a maximum value (cu,) can be 
calculated by: w 

a v (P, S. w) = - a **~ (S) 



9. 



P 1/p (E[li*V' P 

where 

°W(S) is a base value, dependent on the transformation used S and on the (base) sub-band considered for 
the insertion, of the maximum weighting coefficient ensuring imperceptibility when a single coefficient of said 
sub-band is modulated, 
P is strictly greater than 2, and 

giiviP] is the mathematical expectation of the function IvvlP. 

Method of determining a partitioning according to one of Claims 2 or 3, characterised in that the detectability 
criterion is calculated by using, for each coefficient (Xj) to be modulated, a weighting coefficient in accordance with 
a law of the form a s = kj.a v , where kj is a modulation factor dependent on the coefficients situated close to the 
relevant coefficient (X s ) on said region and a v is equal to said maximum value of the weighting coefficient. 

Method of determining a partitioning according to one of claims 1 to 4, characterised in that it comprises a step 
(E307) of applying a distortion before calculating the detectability criterion. 

Method of determining a partitioning according to one of Claims 1 to 5, characterised in that it includes a prior step 
(E105) of comparing the size (P) of each region (R) of the partitioning with a minimum size (P stat ) corresponding 
to the minimum size with a statistically significant sample for the detection of a watermarking signal inserted in 
said region. ~ 

Method of determining a partitioning according to one of Claims 1 to 6, characterised in that, for each region (R) 
of the partitioning, all the coefficients of said region are modulated by inserting the watermarking signal (W) a 
detectability amplitude (T) is calculated from said detectability criterion and said detectability criterion is validated 
by comparing said detectability amplitude (T) with a predetermined threshold value (Tc). 

Method of determining a partitioning according to claim 7, characterised in that said detectability criterion is vali- 
dated if a minimum modulation length (P min (x)) calculated for ensuring detectability of the inserted watermarking 
signal (W) is less than or equal to the cardinal number of said region (R). 

Method of determining a partitioning according to claim 1 , characterised in that it comprises the following steps: 

- generation (E501 ) of a centred pseudo-random sequence (w) equal in size to the cardinal number (N) of a set 
of coefficients representative of at least a part of an image (I), formed from centred pseudo-random sub- 
sequences (w*); 

- modulation (E502) of the said set of coefficients by the said centred pseudo-random sequence (w) in order to 
insert the same information bit on the said set of coefficients; and 

- checking (E506, E507, E514, E516), over each region, of a criterion of detectability of the information bit 
inserted by modulation. 

1 0. Method of determining a partitioning according to Claim 9, characterised in that it comprises a single step of applying 
a distortion (E503) to the said set of modulated coefficients, before the step of checking (E506, E507 E51 4 E51 6) 
of a detectability criterion. 

11. Method of determining a partitioning according to Claim 10, characterised in that the said set of coefficients is a 
set of spectral coefficients (X) originating from a spectral transformation of a digital image, and in that the said step 
(E503) of applying a distortion includes the following sub-steps: 

- inverse spectral transformation (E503a) of the said set of modulated coefficients in order to obtain a water- 
marked image; 

- application (E503b) of a distortion to the said watermarked image; and 
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- spectral transformation (E503c) of the watermarked image in order to obtain a set of modulated coefficients 
after distortion. 

12. Method of determining a partitioning according to one of Claims 9 to 11, characterised in that the centred pseudo- 
random sequence (w) is formed from pseudo-random sub-sequences (wk) of identical size. 

13. Method of determining a partitioning according to one of Claims 9 to 12, characterised in that the pseudo-random 
sub-sequences (wk) have a size greater than or equal to a minimum number of coefficients forming a region of 
the said partition of a size which is statistically significant for validating the said criterion of detectability over the 
said region. 

14. Method of determining a partitioning according to one of Claims 9 to 13, characterised in that, at the generation 
step (E501 ), each centred pseudo-random sub-sequence (wk) is created by generating half of the pseudo-random 
samples of the said sub-sequence via a pseudo-random number generator and by symmetrising the said pseudo- 
random samples generated. 

15. Method of determining a partitioning according to Claim 14, characterised in that, at the generation step (E501 ), 
each centred pseudo-random sub-sequence (wk) is created by moreover permutating the said pseudo-random 
samples using a secret key. 

16. Method of determining a partitioning according to one of Claims 9 to 13, characterised in that, at the generation 
step (E501), each centred pseudo-random sub-sequence (wk) is created by generating pseudo-random random 
samples via a pseudo-random number generator and by redistributing over each sample the sum of the said 
pseudo-random samples generated. 

17. Method of determining a partitioning according to one of Claims 7 to 16, characterised in that a partitioning by 
iterative division (E103) of the digital image is effected and in that, for each region at a rank of said partitioning, a 
partitioning of a directly higher rank of said region is validated if and only if there exist at least two sub-regions of 
said region for which said detectability criterion is validated. 

1 8. Method of determining a partitioning according to Claim 1 7, characterised in that, for each region at a rank of said 
partitioning, the value of a capacity (Q) of the digital image is incremented (E11 2) when the partitioning of directly 
higher rank is not validated, said capacity (Q) being equal to the cardinal number of the distinct regions obtained 
by said partitioning on which the detectability criterion is validated. 

19. Method of determining a partitioning according to Claim 18, characterised in that, for each region at a rank of said 
partitioning, the value of a capacity (Q) is incremented (E112) when the partitioning of immediately high rank is 
not validated, said value of the capacity (Q) is compared (E11 5) with a fixed capacity value (Qs) and the partitioning 
is reiterated for a region having the highest detectability amplitude amongst all the other regions to be processed 
when said capacity value is lower than said fixed capacity value. 

20. Method of determining a partitioning according to one of Claims 7 to 16, characterised in that a partitioning by 
iterative merging of the digital image is effected and in that, for each region of said partitioning for which said 
detectability criterion is not validated, said region is merged with another region of said partitioning. 

21 . Method of determining a partitioning according to Claim 20, characterised in that said other region of the partitioning 
is if possible a region for which said detectability criterion is not validated. 

22. Method of inserting, into a digital image, a watermark comprising one or more information bits, characterised in 
that it comprises a method of determining a partition according to one of Claims 1 to 21, and a step (E520) of 
inserting the information bit or bits by modulation of the coefficients of one or more regions of the said partition. 

23. Method of decoding a watermarking signal in a digital image, inserted by modulating coefficients (X) representing 
said image in distinct regions forming a partitioning of said digital image, characterised in that a partitioning into 
distinct regions of the digital image to be decoded is determined by effecting an adaptive partitioning as a function 
of a criterion of detectability of an information bit obtained by demodulating coefficients on each region (R). 

24. Decoding method according to Claim 23, characterised in that it includes a prior step (E201 ) of comparing the size 
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(P) of each region of the partitioning with a minimum size (P stat ) corresponding to the minimum size of a statistically 
significant sample for the detection of a watermarking signal inserted in said region. 

>. Decoding method according to one of Claims 23 or 24, characterised in that, for each region (R) of the partitioning 
a detection amplitude (T) is calculated (E203) from said detectability criterion and said detectability criterion is 
validated (E203) by comparing said detection amplitude (T) with a predetermined threshold value (Td) for decoding. 

L Decoding method according to Claim 25, characterised in that the predetermined threshold value (Td) for decoding 
is less than said predetermined threshold value (Tc) used during the partitioning determination method accordinq 
to Claim 7. 3 



27. Device for determining a partitioning into distinct regions of a digital image for inserting a watermarking signal by 
modulating coefficients of said regions representing a digital image, characterised in that it comprises partitioning 
means (500, 502, 503) adapted to effect an adaptive partitioning as a function of a criterion of detectability of an 
information bit inserted on each region. 

28. Partitioning determination device according to claim 27, characterised in that it also includes means (500, 502, 
503) for applying a distortion before calculating the detectability criterion. 

29. Partitioning determination device according to one of claims 27 to 28, characterised in that it also has means (500, 
502, 503) of comparing the size of each region of the partitioning with a minimum size corresponding to the minimum 
size of a statistically significant sample for the detection of a watermarking signal inserted in the said region. 

30. Partitioning determination device according to one of Claims 27 to 29, characterised in that the partitioning means 
(500, 502, 503) comprise means (500, 502, 503) of modulating all the coefficients of each region by inserting the 
watermarking signal, means (500, 502, 503) of calculating a detectability amplitude from said detectability criterion 
and means (500, 502, 503) of validating said detectability criterion, adapted to compare said detectability amplitude 
with a predetermined threshold value. 

31. Partitioning determination device according to claim 27, characterised in that it includes: 

- means (500, 502, 503) for generating a centred pseudo-random sequence (w) equal in size to the cardinal 
number (N) of a set of coefficients representative of at least a part of an image (I), formed by centred pseudo- 
random sub-sequences (w k ); 

- means (500, 502, 503) for modulation of the said set of coefficients by the said centred pseudo-random se- 
quence (w) in order to insert the same information bit on the said set of coefficients; and 

- means 500, 502, 503) for checking, over each region, a criterion of detectability of the information bit inserted 
by modulation. 

32. Partitioning determination device according to Claim 31 , characterised in that it comprises means (500, 502, 503) 
for applying a distortion to the said set of modulated coefficients. 

33. Partitioning determination device according to Claim 32, characterised in that the said set of coefficients is a set 
of spectral coefficients (X f ) originating from a spectral transformation of a digital image, and in that the said means 
for applying a distortion comprise: 

- means (500, 502, 503) for inverse spectral transformation of the said set of modulated coefficients in order to 
obtain a watermarked image; 

• means (500. 502, 503) for application of a distortion to the said watermarked image; and 

- means (500, 502, 503) for spectral transformation of the watermarked image in order to obtain a set of mod- 
ulated coefficients after distortion. 

34. Partitioning determination device according to one of Claims 30 to 33, characterised in that the partitioning means 
(500, 502, 503) are adapted to effect a partitioning by iterative division of the digital image and in that, for each 
region at a rank of said partitioning, the validation means are adapted to validate a partitioning of a directly higher 
rank of said region if and only if there exist at least two sub-regions of said regions for which said detectability 
criterion is validated. 
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35. Partitioning determination device according to Claim 34, characterised in that it also comprises means (500, 502, 
503) of incrementing the value of a capacity of the digital image adapted, for each region at a rank of said parti- 
tioning, to increment the said capacity value when the partitioning of directly higher rank is not validated, said 
capacity being equal to the cardinal number of the distinct regions obtained by said partitioning on which the 
detectability criterion is validated. 



36. Partitioning determination device according to Claim 34, characterised in that it also comprises incrementation 
means (500, 502, 503) adapted, for each region at a rank of said partitioning, to increment the value of a capacity 
when the partitioning of directly higher rank is not validated, comparison means (500, 502. 503) adapted to compare 
said value of the capacity with a fixed capacity value and reiteration means (500, 502, 503) adapted to reiterate 
the partitioning for a region having the highest detectability amplitude amongst all the other regions to be processed 
when said capacity value is lower than said fixed capacity value. 

37. Partitioning determination device according to one of Claims 30 to 33, characterised in that the partitioning means 
(500, 502, 503) are adapted to effect a partitioning by iterative merging of the digital image and are adapted, for 
each region of said partitioning for which said detectability criterion is not validated, to merge said region with 
another region of said partitioning. 

38. Partitioning determination device according to Claim 37, characterised in that said other region of the partitioning 
is if possible a region for which said detectability criterion is not validated. 

39. Partitioning determination device according to one of Claims 27 to 38, characterised in that it is incorporated in a 
microprocessor (500), a read only memory (502) being adapted to store a program for determining a partitioning 
and a random access memory (503) comprising registers adapted to store variables modified during the running 
of said program. 



40. Device for inserting, into a digital image, a watermark comprising one or more information bits, characterised in 
that it comprises a device (1 1 ) for determining a partition according to one of Claims 27 to 39, and means (1 2) for 
inserting the information bit or bits by modulation of the coefficients of one or more regions of the said partition. 

41. Insertion device according to Claim 40, characterised in that it is incorporated into a microprocessor (500), a read- 
only memory (502) being suitable for storing in memory a program for inserting a watermark by modulation, and 
a random-access memory (503) including registers suitable for storing in memory the variables modified during 
the running of the said program. 

42. Device for decoding a watermarking signal in a digital image, inserted by modulation of coefficients representing 
said image in distinct regions forming a partitioning of said digital image, characterised in that it comprises means 
(500, 502, 503) of partitioning of the digital image to be decoded into distinct regions, adapted to effect an adaptive 
partitioning as a function of a criterion of detectability of an information bit obtained by demodulation of the coef- 
ficients on each region. 

43. Decoding device according to Claim 42, characterised in that it also has means (500, 502, 503) of comparing the 
size of each region of the partitioning with a minimum size corresponding to the minimum size of a statistical 
significant sample for the detection of a watermarking signal inserted in said region. 

44. Decoding device according to one of Claims 42 or 43, characterised in that it comprises means (500, 502, 503) of 
calculating a detection amplitude from said detectability criterion and means (500, 502, 503) of validating said 
detectability criterion adapted to compare said detection amplitude with a threshold value predetermined for the 
decoding. 

45. Decoding device according to Claim 44, characterised in that the predetermined threshold value for the decoding 
is less than said predetermined threshold value used during the partitioning determination method according to 
Claim 7. 

46. Decoding device according to one of Claims 42 to 45, characterised in that it is incorporated in a microprocessor 
(500), a read only memory (502) being adapted to store a program for decoding a watermarking signal and a 
random access memory (503) comprising registers adapted to store variables modified during the running of said 
program. 
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47. Method of processing an image, characterised in that it comprises the following steps: 

inputting an image (E500); 

inserting a first watermarking information into said inputted image (E502); 
applying a distortion to the image obtained in said insertion step (E503a-c); 
• detecting said first watermarking information from the image obtained' in said distortion applvinq step 
(E504-E508); and a K 

- determining an amount of a second watermarking information which can be inserted into said inputted image 
in accordance with the result of the detection obtained in said detection step (E510). 

48. Processing method according to Claim 47, characterised in that it futher comprises the step of displaying the 
amount of the watermarking information determined in said determining step. 

49. Image processing apparatus, characterised in that it comprises: 

means (500. 502, 503) for inputting an image; 

- means (500, 502, 503) for inserting a first watermarking information into said inputted image; 

- means (500, 502, 503) for applying a distortion to the image obtained by said insertion means; 

- means (500. 502, 503) for detecting said first watermarking information from the image obtained by said dis- 
tortion applying means; and 

- means (500, 502, 503) for determining an amount of a second watermarking information which can be inserted 
into said inputted image in accordance with the result of the detection obtained by said detection means. 

50. Image processing apparatus according to Claim 49, characterised in that it further comprises the means (500, 502, 
503) for displaying the amount of the watermarking information determined by said determining means. 

51. Method of inserting a watermarking signal (w) in a set of coefficients (X) representing a digital image (I), in which 
at least one subset of coefficients is modulated according to a linear model using a weighting coefficient (a) the 
weighting coefficient (a) being greater than a minimum value (oq) determined as a function notably of the length 
(P) of the watermarking signal (w) and a detection probability level (p), characterised in that it includes: 

- a step (E602) of calculating a maximum value (a v (N)) of the weighting coefficients functional to the length (P) 
of the watermarking signal (w) ensuring imperceptibility of the watermarking signal. 

52. Watermarking signal insertion method according to Claim 51, characterised in that it also includes: 

- a step (E603) of determining compatible values of the length (P) of the watermarking signal (w) of the weighting 
coefficient (a) and of th e detection probability level (p) so that the value of the weighting coefficient determined 
(a) is less than or equal to said maximum value (ov) and greater than or equal to said minimum value (ct D ) for 
said given values of the length (P) of the watermarking signal (w) and the detection probability level (p). 

53. Watermarking signal insertion method according to one of Claims 51 or 52, characterised in that, during the step 
(E602) of calculating the said maximum value (a v ), a function is used which is independent of said set of coefficients 
to be watermarked (X). 

54. Watermarking signal insertion method according to one of Claims 51 to 53, characterised in that said coefficients 
are spatio-frequency coefficients obtained by spatio-frequency transformation of an image (I), the calculation of 
said maximum value (eg of the weighting coefficient depending on the type of transformation used (S). 

55. Watermarking signal insertion method according to Claim 54, characterised in that said calculation function used 
is of the type 



0L V {P. S, IV): a ^- (S) 



P 1/P (E[lwl P ]) 1/p 
where 
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°We( s ) is a b as© value dependent on the transformation used (S) and on the (base) sub-band considered 
for insertion, of the maximum weighting coefficients ensuring imperceptibly when a single coefficient of said 
sub-band is modulated, 
p is strictly greater than 2, and 

EflwlP] is the mathematical expectation of the function lu/iP. 

56. Watermarking signal insertion method according to one of Claims 51 to 55, characterised in that said minimum 
value (a D ) is determined by using a calculation function also dependent on the variance (a x 2 ) of the coefficients 
to be modulated, the distribution of the watermarking signal (w) and possibly the variance (o n 2 ) of an additive noise 
decorrelated from the coefficients modelling any distortions applied to said coefficients. 

57. Watermarking signal insertion method according to Claim 56, characterised in that said calculation function is of 
the type 



a D (P) = 



JbP- 



where c x 2 is equal to the variance of the coefficients to be modulated, a n 2 is equal to the variance of the 
additive noise and a, b and c are constants which depend on the distribution of the watermarking signal (w) and 
the required detection probability level (p). 

58. Watermarking signal insertion method according to one of Claims 51 to 57, characterised in that the watermarking 
signal (w) is a pseudo-random sequence with a predetermined distribution with a zero mean. 

59. Watermarking signal insertion method according to one of Claims 51 to 58, characterised in that it includes the 
following steps: 

- calculating (E602) a maximum value (ot^N)) of the weighting coefficients for a length (P) of the watermarking 
signal equal to the cardinal number (N) of the set of modulatable coefficients (X) ensuring imperceptibly of 
the watermarking signal (w) at the limit of perceptibility; 

- calculating (E603) the length (P°) of the watermarking signal (w) for a minimal value (Oq) of the weighting 
coefficient equal to said maximum value calculated (<x v (N)) and for a predetermined detection probability level 
(p); and 

- comparing (E606) said calculated length (P°) with the cardinal number (N) of the set of modulatable coefficients 
(X). 

60. Watermarking signal insertion method according to Claim 59, characterised in that, if at the comparison step (E606) 
said calculated length (P°) is less than or equal to the cardinal number (N) of the set of modulatable coefficients 
(X), there is inserted a watermarking signal (w) of length (P) at least equal to said calculated length (P°) by mod- 
ulating a subset of coefficients with a cardinal number equal to said length (P) according to a linear model using 
a weighting coefficient (a) equal to said calculated maximum value (a v (N)). 

61 . Watermarking signal insertion method according to Claim 59, characterised in that, if at the comparison step (E606) 
the said length (P°) is greater than the cardinal number (N) of the set of modulatable coefficients (X), there is 
calculated the detection probability level (p') obtained for a minimum value (cc D ) of the weighting coefficient equal 
to said calculated maximum value (a v (N)) for a length (P) of the watermarking signal (w) equal to the cardinal 
number (N) of the set of modulatable coefficients (X). 

62. Watermarking signal insertion method according to one of Claims 51 to 58, characterised in that it includes a step 
of calculating a threshold value (P*) of the length of the watermarking signal (w) determined so that said minimum 
value (ceo) of the weighting coefficients is equal to said maximum value (a v ) of the weighting coefficient for a length 
(P) of the watermarking signal (w) equal to said threshold value (P*) and a predetermined detection probability 
level (p). 

63. Device for inserting a watermarking signal (w) in a set of coefficients representing a digital image (I), in which at 
least one subset of coefficients is modulated according to a linear model using a weighting coefficient (a), the 
weighting coefficient being greater than a minimum value (oq) determined as a function notably of the length (P) 
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ol the watermarking signal (w) and a detection probability level (p), characterised in that it comprises means (500 
502, 503) of calculating a maximum value (<x v ) of the weighting coefficient as a function of the length (P) of the 
watermarking signal ensuring imperceptibiiity of the watermarking signal (w). 

64. Watermarking signal insertion device according to Claim 63, characterised in that it also comprises means (500 
502, 503) of determining compatible values of the length (P) of the watermarking signal, of the weighting coefficient 
(a) and of the detection probability level (p) so that the value of the weighting coefficient determined (a) is less 
than or equal to said maximum value (a v ) and greater than or equal to said minimum value (oto) for said determined 
values of the length (P) of the watermarking signal and of the detection probability level (p). 

65. Watermarking signal insertion device according to one of Claims 63 or 64, characterised in that the means (500 
502, 503) of calculating said maximum value (a v ) use a function independent of said set of coefficients to be 
watermarked. 

66. Watermarking signal insertion device according to one of Claims 63 to 65, characterised in that it comprises means 
(500, 502, 503) for the spatio-frequency transformation of an image so that the coefficients are spatio-frequency 
coefficients obtained by spatio-frequency transformation of an image, the calculation of said maximum value (a v ) 
of the weighting coefficient depending on the type of transformation used (S). 

67. Watermarking signal insertion device according to Claim 66, characterised in that said calculation means (500 
502, 503) use a function of the type 

<M P - S ' *>= l/B 0 ^? 1/B 

P 1/P (E[lwlV /p 

where 

a base(S) is a base value dependent on the transformation used S and on the (base) sub-band considered for 
the insertion, of the maximum weighting coefficient ensuring perceptibility when a single coefficient of the said 
sub-band is modulated, 
P is strictly greater than 2. and 

E[ltvlP] is the mathematical expectation of the function IwlP. 

68. Watermarking signal insertion device according to one of Claims 63 to 67, characterised in that said minimum 
value (a D ) is determined by using a calculation function also depending on the variance (o/) of the coefficients 
to be modulated, the distribution of the watermarking signal (w) and possibly the variance (c n 2) of an additive noise 
decollated from the coefficients modelling any distortions applied to said coefficients. 

69. Watermarking signal insertion device according to Claim 68, characterised in that said calculation function is of 
the type 



where o x 2 is equal to the variance of the coefficients to be modulated, o n 2 is equal to the variance of the 
additive noise and a, b and c are constants which depend on the distribution of the watermarking signal and the 
required detection probability level. 

70. Watermarking signal insertion device according to one of Claims 63 to 69, characterised in that the watermarking 
signal (w) is a pseudo-random sequence with a predetermined distribution with a zero mean. 

71 . Watermarking signal insertion device according to one of Claims 63 to 70, characterised in that it also comprises: 

- means (500, 502, 503) of calculating the length (P°) of the watermarking signal for a minimum value (Oq) of 
the weighting coefficient equal to a maximum value (a v (N)) calculated for a length (P) of the watermarking 
signal equal to the cardinal number (N) of the set of modulatable coefficients ensuring imperceptibiiity of the 
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watermarking signal at the limit of perceptibility and for a predetermined detection probability level (p); and 
- means (500, 502, 503) of comparing said calculated length (P°) with the ca.dinai number (N) of the set of 
modulatable coefficients (X). 

72. Watermarking signal insertion device according to Claim 71, characterised in that it also comprises means (500, 
502, 503) of inserting a watermarking signal (w) with a length (P) at least equal to said calculated length (P°), 
adapted to modulate a subset of coefficients with a cardinal number equal to said length (P) according to a linear 
model using a weighting coefficient (a) equal to said calculated maximum value (a v (N), if said calculated length 
(P°) is less than or equal to the cardinal number (N) of the set of modulatable coefficients (X). 

73. Watermarking signal insertion device according to Claim 71, characterised in that it also comprises means (500, 
502, 503) of calculating a detection probability level (p') obtained for a minimum value (a D ) of the weighting coef- 
ficient equal to said maximum value (a v (N)) calculated for a length (P) of the watermarking signal equal to the 
cardinal number (N) of the set of modulatable coefficients. 

74. Watermarking signal insertion device according to one of Claims 63 to 70, characterised in that it comprises means 
(500, 502, 503) of calculating a threshold value (P*) of the length of the watermarking signal determined so that 
said minimum value (<*□) of the weighting coefficient is equal to said maximum value (a v ) of the weighting coefficient 
for a length (P) of the watermarking signal (w) equal to said threshold value (P*) and a predetermined detection 
probability level (p). 

75. Watermarking signal insertion device according to one of Claims 63 to 74, characterised in that it is incorporated 
in a microprocessor (500), a read only memory (502) storing a program for inserting a watermarking signal in an 
image and a random access memory (503) containing registers adapted to store variables modified during the 
running of said program. 

76. Computer, characterised in that it comprises means adapted to implement a partitioning determination method 
according to one of Claims 1 to 21 . 

77. Computer, characterised in that it comprises means adapted to implement the method of inserting a watermark 
according to Claim 22. 

78. Computer, characterised in that it comprises means adapted to implement the decoding method according to one 
of Claims 23 to 26. 

79. Computer, characterised in that it comprises means adapted to implement a method of processing an image ac- 
cording to one of Claims 47 to 48. 

80. Computer, characterised in that it comprises means adapted to implement a method of inserting a watermarking 
signal according to one of Claims 51 to 62. 

81. Digital image processing apparatus, characterised in that it comprises means adapted to implement a partitioning 
determination method according to one of Claims 1 to 21 . 

82. Digital image processing apparatus, characterised in that it comprises means adapted to implement the method 
of inserting a watermark according to Claims 22. 

83. Digital image processing apparatus, characterised in that it comprises means adapted to implement the decoding 
method according to one of Claims 23 to 26. 

84. Digital image processing apparatus, characterised in that it comprises means adapted to implement a method of 
processing an image according to one of Claims 47 or 48. 

85. Digital image processing apparatus, characterised in that it comprises means adapted to implement a method of 
inserting a watermarking signal according to one of Claims 51 to 62. 

86. Digital printer, characterised in that it comprises means adapted to implement a partitioning determination method 
according to one of Claims 1 to 21 . 
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87. Digital printer, characterised in that it comprises means adapted to implement the method of inserting a watermark 
according to Claim 22. 

88. Digital printer, characterised in that it comprises means adapted to implement the decoding method according to 
one of Claims 23 to 26. 

89. Digital printer, characterised in that it comprises means adapted to implement a method of processing an image 
according to one of Claims 47 or 48. 

90. Digital printer, characterised in that it comprises means adapted to implement a method of inserting a watermarking 
signal according to one of Claims 51 to 62. 

91. Digital photographic apparatus, characterised in that it comprises means adapted to implement a partitioning de- 
termination method according to one of Claims 1 to 21. 

92. Digital photographic apparatus, characterised in that it comprises means adapted to implement the method of 
inserting a watermark according to Claim 22. 

93. Digital photographic apparatus, characterised in that it comprises means adapted to implement the decoding meth- 
od according to one of Claims 23 to 26. 

94. Digital photographic apparatus, characterised in that it comprises means adapted to implement a method of 
processing an image according to one of Claims 47 or 48. 

95. Digital photographic apparatus, characterised in that it comprises means adapted to implement a method of in- 
serting a watermarking signal according to one of Claims 51 to 62. 

96. Digital camera, characterised in that it comprises means adapted to implement a partitioning determination method 
according to one of Claims 1 to 21 . 

97. Digital camera, characterised in that it comprises means adapted to implement the method of inserting a watermark 
according to Claim 22. 

98. Digital camera, characterised in that it comprises means adapted to implement the decoding method according to 
one of Claims 23 to 26. 

99. Digital camera, characterised in that it comprises means adapted to implement a method of processing an image 
according to one of Claims 47 or 48. 

100. Digital camera, characterised in that it comprises means adapted to implement a method of inserting a watermark- 
ing signal according to one of Claims 51 to 62. 

101 .Computer program product, characterised in that it is adapted to implement a partitioning determination method 
according to one of Claims 1 to 21 . 

102. Computer program product, characterised in that it is adapted to implement a method of inserting a watermark 
according to Claim 22. 

103. Computer program product, characterised in that it is adapted to implement a decoding method according to one 
of Claims 23 to 26. 

104. Computer program product, characterised in that it is adapted to implement a method of processing an image 
according to one of Claims 47 or 48. 

105. Computer program product, characterised in that it is adapted to implement a method of inserting a watermarking 
signal according to one of Claims 51 to 62. 

106.lnformation storage means readable by a microprocessor, characterised in that it stores a program implementing 
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a partitioning determination method according to one of Claims 1 to 21 . 

107. Information storage means readabie by a microprocessor, characterised in that it stores a program implementing 
a method of inserting a watermark according to Claim 22. 

108.lnformation storage means readable by a microprocessor, characterised in that it stores a program implementing 
a decoding method according to one of Claims 23 to 26. 

1 (^.Information storage means readable by a microprocessor, characterised in that it stores a program implementing 
a method of processing an image according to one of Claims 47 or 48. 

HO.Information storage means readable by a microprocessor, characterised in that it stores a program implementing 
a method of inserting a watermarking signal according to one of Claims 51 to 62. 
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